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APPLICATIONS OF TWO OSCULATORY FORMULAS 


By Joun L. RoBERTS 


INTRODUCTION 


The main purpose of this paper is to illustrate how Mr. Jenkins’ osculatory 
formulas’ (A) and (B) can be applied in a convenient manner. The first section 
of this paper will be little more than a summary of some of the formulas con- 


tained in the other three articles. The second section will contain the appli- 
cations. 


I. SOME MATHEMATICS OF THE FORMULAS 


The Woolhouse notation will in this paper be used to stand for the differences 
of uz, Which represents the given values of afunction. The general formulas are 


Ye = Yo + cAyo + 3a(x — 1)B + ga(x — I) — 3)C; (1) 
and 
Yr = Uy + xa, + 3a(a — 1)B + $a(a — 1)(a — 3)C. (2) 
The special formulas belonging to (2) are 
B=b-— jdandC = q — je, (A) 
where 6 and d are defined by b = 4(bo + bi) and by d = 3(d) + di); and 
B= bandC = 0. (B) 


The special formulas belonging to (1) are 


Yo = Uo + tbo, b, and 5 0; (C) 


and 
Yo = Uy — sedo ; B=b6- dd, and C a > bey ‘ (D) 


Formula (C) is equivalent to Mr. Jenkins’ formula (A). Also (D) is equivalent 
to his formula (B). 


1 This paper presupposes a knowledge of three other articles. The first one by Mr. 
Wilmer A. Jenkins is entitled ‘‘Graduation Based on a Modification of Osculatory Inter- 
polation,’ and is printed in the October 1927 issue of the Transactions of the Actuarial 
Society of America. The other two papers are mine. One of them is entitled ‘Some 
Practical Interpolation Formulas,” and is printed in the September 1935 issue of these 
Annals. The other one entitled ‘‘A Family of Osculatory Formulas”’ is printed in the 
October 1935 issue of the Transactions. 
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II. APPLICATIONS OF (C) AND (D) 


First, there is the problem of selecting suitable examples to which (C) and (D) 
can be applied. Secondly, we will then apply in a convenient manner the 
formulas to these examples. 

The problem of selecting suitable examples will now be considered. ‘‘The 
non-reproducing characteristic of’? formula (D) ‘“‘raises the question of what 
will happen in the graduation of a series whose fourth differences are all posi- 
tive, say. The answer is that the graduated series will lie everywhere below 
the observed points and that the observations will not be correctly represented 
by the interpolated series.”” On the other hand, if we select a series whose 
fourth differences change frequently in sign, (D) because of its non-reproducing 
characteristic has valuable smoothing possibilities. In like manner, (C) may 
be valuable when the second differences change frequently in sign. Mr. Jenkins 
gives at quinquennial ages rates of mortality which were graphically determined 
from the published American Men Ultimate Experience. Since the fourth 
differences of these rates change frequently in sign, we will apply (D) to a few 
of these rates. So far as I know no suitable actuarial examples have been 
found to which (C) can be applied. However, there is the possibility that (C) 
might be valuable in some sciences. Since I do not know of any suitable real 
example to which (C) can be applied, we will apply it to a trivial series whose 
second differences change frequently in sign. 

We are now ready to apply in a convenient manner (C) and (D) to the 
examples selected in the preceding paragraph. 

First, we will apply (C). I have in my other article applied (B) in a con- 
venient manner. This method with little change can be applied to (C). If 
it is desired to apply (C) at either end of the table where values of uz are not 
available for the calculation of the second differences, it can be assumed they 
vanish. It is convenient if S and S’ represent respectively the major differ- 
ences Au, and Au; in such a manner that they are arranged centrally in the 
working illustration. It is convenient if s and s’ represent respectively the 
minor differences dy, and dy,. The quantity yo can be computed by yo = 
Uo + 4bo , and y; can be computed in like manner. Since we wish in the working 
illustration of (C) to interpolate four values between yo and y, the middle 
s = dy4 = .2Ay, ands’ = .04B = .02(b) + b:). We can by the use of the 
foregoing method apply (C) to suitable functions, whose given values can be 
represented by f(r). Then, it follows from the definition of uz that f(r) = uz. 
It might prevent confusion if it is stated that x and r are related to each other 
in such a way that we always interpolate between yo and y,;. We shall now 
apply (C) to the case when f(r) represents the trivial series shown at top of 
page 3. 

Finally, we will apply (D). Mr. Henderson has applied (A) in a very con- 
venient manner. His method with little change can be applied to (D). If it 
is desired to apply (D) at either end of the table where values of u, are not 
available for the calculation of the differences required, it can be assumed 
that the fourth differences that can not be computed vanish, and the required 
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f(r) S S? Ye s s? 
— itl heemecaiead _ a 

0 1 0 1.0 | 18 

1 | 118 | 14 

2 | 1 1.32 10 — .04 

3 1.42 06 

4 | 1.48 02 

5 | 2 _2 1.5 | 

6 | 1.5 | 

7 ml 1.5 | 00 00 

~ 1.5 

9 1.5 

10 1 2 1.5 02 

11 1.52 | .06 

12 1 1.58 | 10 04 

13 1.63 14 

14 1.82 18 

15 2 0 2.0 | 

















differences can be filled in consistently with that assumption. 


It is convenient 


° 2 3 . . . 2 
if S, S°, and S° represent respectively the major differences Au,, A’u,, and 
3 * ‘ . ° 
A’u; in such a manner that they are arranged centrally in the working illustra- 
° ° ° ° 2 3 ° ° 
- tion. It is convenient if s, s, and s° represent the minor differences so that 


by definition s = s, = dyz, 8 = sf = Syz-2, ands’ = S yz. 


The first 


s = Sy_2 = .04(bo — 3d). The last s’ = dy. = .04(b, — 4d,). The quan- 
tity yo can be computed by yo = Uo — gedo, and y, can be computed in like 


manner. 


apply (D) to the quinquennial rates of mortality. 











The middle s = 6y4 = .2Ay — s*. We are now in position to 





Age Rate S S? S3 s 

72 .07010 
.03808 

77 . 10818 .00861 
.04669 .01799 

82 . 15487 .02660 — .03745 
.07329 — .01946 

87 . 22816 .00714 . 14518 
.08043 . 12572 

92 .30859 . 13286 .00000 
. 21329 12572 

97 .92188 . 25858 .00000 








Yr 


.15591 


. 168522 
. 182049 
. 196092 
. 210252 


.22413 
. 237590 
. 251567 
. 267260 
. 285868 


.30859 

. 336596 
.370922 
.412574 
.462558 


.92188 
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$s 


12612 
13527 
.014043 
14160 
13878 


13460 
13977 
015693 
18608 
22722 


28006 
34326 
.041652 
49984 
59322 


s? 


.001314 
915 

516 

117 

— .000282 


— .000682 
.000517 
1716 
2915 
4114 


.005314 
6320 
7326 
8332 
9338 


.010343 





.001199 


.001006 





SOME SIMPLE DEVELOPMENTS IN THE USE OF THE 
COEFFICIENT OF STABILITY 


By C. H. ForsytTu 
Some time ago the writer proposed’ a coefficient of stability C, to be used 


to measure the stability of a statistical series, where that coefficient is defined 
by the relation 


iS] 


C. = (1) 


=I 


where M denotes the arithmetic mean and o° the square of the dispersion of 
the terms of the series. It was proposed to regard series as unstable (Lexian) 
for which the value of the coefficient exceeded unity, and stable otherwise. 
The only essential way in which such a procedure differs in results from the 
traditional method is that it includes as stable those series for which the value 
of the coefficient lies between unity and q the probability of failure of the event 
under investigation—series which would be classed as unstable according to 
the traditional method. Stable series—according to either standard—are found 
so rarely in practice and therefore so many series are accepted as fairly stable 
which come anywhere near meeting the requirements that replacing q by unity 
as the line of demarcation affects the classification of no known series but 
adds to the effectiveness of the avowed purpose and use of the proposed coeffi- 
cient—to avoid the round-about work of computing values of probabilities. 
Another merit of the use of the coefficient is that it enables one to measure 
and therefore compare the stability of several series—a feature which we shall 
illustrate later. 

In brief, such a coefficient provides a means of introducing the whole Lexian 
theory into Federal publications such as those on vital statistics, since a com- 
parison of the values of the coefficient for, say different communities or countries, 
would be readily grasped by any reader, whereas the traditional method would 
prove too subtle and laborious, and allow no ready comparison of results. 

For purpose of orientation let us illustrate the situation by analyzing a simple 
series both ways—the traditional way and by the use of the coefficient of sta- 
bility. As an example, let us consider the death rates of white infants under 
one year of age for 1919 (considered on page 89 of the Handbook) for those 
states whose frequencies of births are comparable or which vary little from 


1 Journal of the American Statistical Association, June, 1932. 
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their average of 47,830—where the number of deaths for each state has been 


adjusted to this average as a base. 


Adjusted 
Deaths X 
3350+ 
4700 
3732 
3253 
3686 
ak Sig ch hic atala! dary 3159 
N. Car 3541 
’ 3732 
nee 3780 


FORSYTH 


X — 3659 


-309 
1041 
73 
—406 
27 
-500 
-118 
73 
121 


(X — 3659)? 


95481 
1083681 
5329 
164836 
729 
250000 
13924 
5329 
14641 


9)32933 
M = 3659 


)1633950 
181550 


o 


1335-1333 


The traditional method would be: 
The mean M = np = 3659 where n = 47,830. 


3659 ~ 44171 


He e = = 
ence P = F7g39 8Nd 4 = F730 


and o = npq = 3659 (a) = 3878 


whence og = 58.15 


which is the value of the dispersion we should expect if the basic probability 
were constant throughout. But the value of the dispersion proves to be 
o = 0181550 = 426, and the comparison of the values shows that the basic 
probability to be very variable and therefore the series to be very unstable or 
Lexian. 

The computation of the value of the coefficient of stability is much more 
simple and direct 


C o 181550 


:= > Foe Ci == 49.6 
M 3659 


whose excess over unity also clearly indicates the instability of the series. 
Since proposing the coefficient of stability the writer has been impressed by 
the overwhelming proportion of existing series (such as birth rates, various kinds 
of death rates, etc.) which employ arbitrary bases (such as “per thousand,” 
“per ten thousand,” etc.) usually without mention of the actual base. It is 
obvious, of course, that such rates, or occurrences per arbitrary base, say b, 
can first be adjusted to give occurrences per actual base, say B (assuming that 
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base B* can be determined) but the work can evidently be performed much 
easier. For, since the original series (per arbitrary base b) X,, Xe2,--- Xw 
. B r r , 
would become, on adjustment, b Ais = Xe, see oXy, the mean would become 
B : B \ 
7M and the square of the dispersion 37)? whence the formula for the coeffi- 
b 


cient of stability would become 


o B 
=iF (2) 


As an example, let us consider the general death rates, per 10,000, of New 
Zealand for the years 1921-30. 


C, 


X X — 86 (X — 86)? 
1921 87 
1922 88 
1923 90 
1924 83 
1925 83 
1926 87 
1927 85 
1928 85 
1929 88 
1930 86 
10)862 

M = 86.2 


1 
4 
16 
9 
9 
1 
1 
1 
4 
0 


)46 


| 


6 


This example illustrates the danger of using the coefficient of stability unless 
the series consists of actual occurrences or unless the actual base is given due 
consideration. Without due consideration of the actual base (here the popula- 
tion of New Zealand) one might easily fall into the error of regarding the value 
of the coefficient of stability as 4.6/86.2 and, therefore, the series as very 
stable. But the population of New Zealand is about a million and a half and, 
therefore the true value of the coefficient of stability is 


c. — 46 1,500,000 _ 4 
86.2 10,000 

* Strictly speaking, this actual base B should be constant throughout the series; other- 
wise the successive numbers of occurrences—the terms of the series—would not be com- 
parable. Where, however, the base B varies little from term to term—as usually happens 
even in the best of series, such as a series of some kind of rates of the same community 
over a short interval—the variation can be ignored, in which case base B (to which the 
terms of the series are adjusted) usually means the arithmetic mean of the different bases. 
In the first treated above, the investigation was limited to certain states in an effort to 
comply with the rule just mentioned but the example is a poor one since the variations 
are still dangerously too large. The situation is saved by the conclusive results. 











8 Cc. H. FORSYTH 
which shows the series to be unstable. However, before we condemn New 
Zealand’s death rates too severely, let us compare her record with those of 


other important countries, including our own, for the same period. 


General Death Rates (per 10,000) 


M Cc, 
Ne ea ia a eG ke whe Ka wwe Ran 86.2 8 
silk Sako RO ce ch ci dice Se ee ie eae 94.3 90 
I Sh tie ae ee Saeed ane adh hay mn a his ee 120.4 96 
EE bid a ge ee a ek ERE REE RES 137.3 139 
TE llc eink Sidhe cle ee ie he waa kes 151.1 536 
LOCI EE RE, 830 
I Sicko olin OF ced Wah eceeN A EERSS 121.3 1117 
TE a le a ag 2 ea hea aA Re ee 170.3 1129 
Rit ees oan bares len beedses Yaka aad 193.7 2190 
SN eats Paceline ch fans bak B va Goh tg x a Ah cl > BOHR MLA 163.5 2760 
EE ee ES Tee ee ee ee 125.4 6040 
Ae Oeil is hi Ses tired hah en ia ANN MG 206.4 6800 


These results show how extremely unstable most series of general death 
rates are and that the series for New Zealand, while unstable according to 
our strict criterion, enjoys quite an enviable position practically in a class by 
itself. Parenthetically, these results also illustrate fairly well the triviality, 
with respect to results, of replacing q by unity as the critical value of the coeffi- 
cient of stability, discussed at the beginning of this article. 

The values of the coefficient listed above would, of course, be reduced some- 
what in most cases if the trend of the series were first eliminated but the writer 
has gone though all this work and found it not worth while—that is, the series 
would still remain markedly unstable. 

Another development proves useful when, as frequently happens, the actual 
base B is unknown to a degree of accuracy desirable for use in formula (2). 


2 
Se «i 


From the inequality I 


we obtain 


Bs (3) 


which is to be used to show how small an actual base should be for the given 
series to be stable. As an example, let us consider the maternal mortality, 
per 10,000 live births, in the so-called expanding registration area of the United 
States. 
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Maternal Deaths in the United States (per 10,000 live births) (Expanding 
Registration Area) 


x X — 66 (X — 66)? 
1923 67 1 1 
1924 66 0 0 
1925 65 —1 1 
1926 66 0 0 
1927 65 -1 1 
1928 69 3 9 
1929 70 4 16 
1930 67 1 1 
1931 66 0 0 
1932 64 —2 4 
10)665 9—4 )33 

66.5 3.3 


Hence, by formula (3), B < o (10,000) or about 200,000. The number 
of live births varies so greatly that we should probably find it impossible to 
agree upon a satisfactory number’ to use as an actual base for such an “ex- 
panding area” but we should all agree that it would be so much greater than 
200,000 that the instability of the series would be unquestioned. 

One must be careful in comparing the results of two or more investigations 
like the one just conducted. For example, the analogous result for Canada 
for the same period yields B < 113,000 and we might conclude, too hastily, 
that the United States series is more stable (or less unstable) whereas any 
knowledge whatever of the numbers of live births of the two countries would 
show that Canada comes much closer to fulfilling her requirement than the 
United States and that the palm must go to Canada. For one thing, Canada 
has about the population of New York city and New York city has about 
100,000 live births annually. In any case, close decisions in matters of this 
kind would be difficult without sufficient information in regard to actual bases. 

There is still another situation which is interesting but of much less impor- 
tance because of the rarity of its occurrence. It will be recalled that the coeffi- 
cient of stability was devised mainly to avoid the use and computation of 
probabilities and that the only difference between the results by the traditional 
method and by the use of the coefficient of stability lies in the trivial replace- 
ment of the critical value qg by unity. In the traditional method of analysis, 
but by comparing the value of the coefficient of stability with q, the coefficient 
is evidently always, strictly speaking, a function of the actual base B. In 
other words, there is no statistical series, however stable it may seem—except 


2 It was in the neighborhood of two million in 1932. 
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for the trivial case when all the terms of the series are exactly the same—but 
what would be unstable if the base were small enough. It is possible to formu- 
late the limit once for all below which the given (otherwise seemingly stable) 
series would prove unstable. 
If, in the relation o” < npg (for stability) we replace p by M/n, g by 1 — M/n 
and then n by B, we obtain 
2 M M 


Ye og cae 
o <M BO B 






















<M-<¢ 











whence, finally 


mM’ 
~ M—o (4) 


V 





B 





where the transference of the term M — o° from one side to the other should 
cause no apprehension since, by hypothesis, « < M and M — o’ is therefore 
always positive. We propose to employ formula (4) in those rare cases where 
the value of the coefficient of stability of actual occurrences—but without 
reference to an actual base—is less than unity—that is, where the given series 
proves to be stable according to the method proposed by the writer—and 
determine the upper limit of the values of the base B for which the series would 
be unstable according to the traditional method of analysis. As an illustra- 
tion, let us consider the familiar series of annual football fatalities in this country 
for the period 1906-1930* (omitting the years when no records were kept). 


Footbali Fatalities 


1906 11 1917 12 
1907 11 1921 12 















1908 13 1923 18 
1909 12 1925 20 
1911 11 1926 9 
1912 13 1927 ; 17 
1913 5 1928 18 
1914 13 1929 12 
1915 15 1930 13 


11.942 
13.055 
Applying formula (4) 


It is easily verified that C, = which is clearly less than unity; whence 


the series clearly seems stable. 












13.055" “ 

B 2 => or ——>D_SCiéCo?': 11588 
13.055 — 11.942 
which shows that the given series is stable as long as the total number of foot- 
ball players exceeds the number 153. A recent news item quoted an estimate 
of the number players participating in games of four hundred colleges as about 
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13,000 and over 600,000 including high schools and all. We can then definitely 
say that the series just considered is stable. Such a conclusion has no bearing, 
of course, upon what might happen if other terms were added to the series. 
It happens that adding the records for the next five years—1931(33), 1932(32), 
1933(27), 1934(25), 1935(30)—would change the whole series to an unstable 
one with C, = 56.9/16.6 = 3.4; but, obviously, the additional records belong 
to a new regime of collection. 














INTERNAL AND EXTERNAL MEANS ARISING FROM THE SCALING 
OF FREQUENCY FUNCTIONS 


By Epwarp L. Dopp 





The scaling’ of frequency functions has been discussed from the standpoint 
of maximum likelihood. But the likelihood criterion to be satisfied sometimes 
leads to a minimum likelihood; and sometimes to neither a maximum nor a 
minimum. Scaling will be studied in this paper with reference to the likelihood 
actually secured, and also with reference to the character of means obtained, 
whether internal or external. 





SECTION 1. INTRODUCTION 


It is well known that a scale obtained in a curve-fitting process is sometimes 
a mean. Thus, with the normal function 













(1) 1 telat 
OO — ’ 
a V 21 
if the scale a is to be obtained from measurements, 71, 22, --- , Xn, We COM- 


monly accept the value 


(2) a= (3 >» wp | ‘ 


that is, the root-mean square of the measurements. Here, the positive value 
of a is naturally taken. It is called the standard deviation, and thought of as 
an appropriate new unit of measure. 

But even with the z’s all negative, and the a taken positive, O. Chisini’ con- 
sidered it proper to regard a as a mean of the 2’s, albeit an external mean. 
From Chisini’s viewpoint, this a whether regarded as positive or negative is 
primarily a solution of 


(3) 


In this sum of squares, the single number a may be substituted for each of the 
x’s. Perhaps this kind of mean should be called a substitutive mean to dis- 
tinguish it from the means of general analysis which are always internal. 















2 2 2 2 2 2 
tate. +a, =~art+at-:--t+a. 









1 Fisher, R. A., “On the mathematical foundation of theoretical statistics,’’ Philo- 
sophical Transactions of the Royal Society of London, Series A, Vol. 222, 309-368, (1921). 
See p. 338. 

2 Chisini, O., ‘‘Sul concetto di media,’’ Periodico di matematico, Series 4, Vol. 9, 106-116, 
(1929). 
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SCALING OF FREQUENCY FUNCTIONS 


The normal function is a particular case of a more general function: 
(4) Constant-a™ e*, o(t) = —t’/p, ¢ = z/a. 


The likelihood method to find the scale a for this function leads to power means, 
including the arithmetic mean, the root-mean-square, root-mean-cube, etc., for 
p = 1, 2, 3, ete. 

The word scale will be used only for a positive number,—which then may be 
regarded as a unit of measurement. 


For measurements, 21, 22, --- , 2, Chisini regarded M as a mean, rclative 
to a function G, provided 


(5) G(x, t2,--- ,%n) = G(M, M,.---, M). 
If a solution of this equation is 

(6) M = F(x, 2, --- , Zn); 

and c is a possible value for the 2’s, it follows at once that 
(7) F(c, c,---,c) =, 


or at least one value of this F is c. Conversely, if (7) is satisfied, it is but a 
change of notation to replace c in (7) by M, and to combine this with (6) to 
obtain 


(8) F(a, 12, --+,%,) = F(M,M, ---, M). 


Hence, this F which in (6) gives explicit form to the implicit M found in (5) 
may also be thought of as a mean-forming function, such as G in (5). Briefly, 
F is a particular G. Thus F(a, x2, +--+ 2n) is a mean of 4%, %2,---,%n, if F 
is so constructed that (7) is satisfied when the arguments are all equal. 
Inasmuch as a frequency function f(t) is non-negative, log, f(t) is real,—say 
¢(t) plus constant. Following R. A. Fisher, it will be convenient to write 
(9) f(t) = Ca e*”, C = Constant 


With location m already determined, the z’s will be thought of as measured. 
from m. And we set 


(10) t= v/a, t; = x;/a, i= as 2, coe, Re 
The “productive” probability—to yield 2, 22, ---, 2,—is then 


(11) L = Iif(t) = C’a * e**". 


This is proportional’ to the “likelihood” of a. Also—it may be noted in 
passing—the productive probability is also proportional to the a posteriori 
probability, if a constant @ priori probability is postulated. The likelihood 
will here be taken as IIf(¢;) itself; and it will be designated by L,—in Fisher’s 


3 Loc. Cit., Fisher, p. 310. 
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notation, L = log Tl. Of course, II and log II take maximum values simul- 
taneously, if at all. From (11) it follows* that 













(12) —a-d log L/aa = n + Zig’ (ti) = Tilt’ (ts) + 1}. 
The equation 
(13) Ltd’ (ti) +n = 0 (i = 1, 2, ---,n) 


will be called the likelihood condition, whether this leads to maximum likeli- 
hood, to minimum likelihood, or to neither. A second differention”’ leads to 


(14) a’-& log L/aa” = Xt’ (ti) — n = DI tid’ (ti) — 1}. 


When negative, this indicates a maximum likelihood; when positive, a minimum 
likelihood for the a obtained from (13). 

Preparatory to the theorems of the next section, just one more matter will 
be discussed. The unit for ¢ is arbitrary; and it may be convenient to write, 
with k ¥ 0, 

(15) 
Then 


(16) 


o(t) = o(ku) = P(u), 










to'(t) = ub'(u). 


Suppose, now, that a positive constant k can be found such that k¢’(k) = —1. 
Then, with ¢ = ku, as postulated, 


(17) 


Thus #’(1) = —1,—or as it will now be written ¢’(1) = —1,—is no more 
restrictive than the condition that some positive k exists such that k¢’(k) = —1. 





1-6’(1) = k@’(k) = —1. 









SECTION 2. GENERAL THEOREMS CONCERNING THE SCALE AS A MEAN 


THEOREM [| 











Given the frequency function 


(18) f(t) = Ca” e*”, é = z/a, t; = z;/a, C = Constant. 
And suppose that | 

(19) ¢’(1) = -1. 

Suppose, also, that for given 2, %2,---, 2,, the likelihood condition (13), 
now written 

(20) Dut (x:/a)¢’(xi/a) + n = 0, 


4 Loc. Cit., Fisher, p. 338. 
5 Loc. Cit., Fisher, p. 339. 
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has a positive solution. 
(21) > F(x, 2, --- » Zn). 


Then this a, the scale, is a mean. 

Proof. With each x; = 0, (20) cannot be satisfied. 

But if, with c # 0, we take each z; = c, and at the same time set a = c, 
then, by (19), 2 = —n; and thus (20) which gives a implicitly is satisfied. 
The explicit a in (21) is therefore such a function F that (7) is satisfied. Hence, 
the scale a is a mean. 


THEOREM II 
Given the frequency function 
(18) f® = Ca* &®, t = x/a, t; = x;/a, C = Constant. 
Suppose that 
(19) (1) = — 1, 
and that 
(22) | tp’(t) | <1 if |¢| <1. 


Moreover, suppose that the likelihood condition (20) for measurements 
21, X2, +++, Xn, has a positive solution a. Then 


(23) a S Maximum |; |. 

Or, suppose that, in place of (22), we have 

(24) l@’@)|>1 if jt] >]; 
and that i¢’(t) keeps the same sign, if |t| > 1. Then 
(25) Minimum | z;| S a. 


Proof. Suppose, if possible, that a > Max |2;|. Then each | 2,;/a| < 1, 
and by (22), | (x:/a)¢’(x:/a) | <1. Then (20) is not satisfied, since | =| < n. 
Thus the hypothesis is contradicted. 

Now (25) is satisfied at once if any x; = 0. But suppose, on the other hand, 
that Min |z;| > 0; and, if possible, that a < Min |2;|. Then, by (24) et 
seq., since | z;/a| > 1, it follows that | =| > n. And thus (20) is again con- 
tradicted. 


THEOREM III 


Given the frequency function 


(18) f@) = Ca" e*, = z/a, t; = z;/a, C = Constant; 
and set y(t) = t¢’(t) + 1. Suppose that 


(26) lim y(t) = a, lim y(t) = B, aB < 0. 


t—0 
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And suppose that ¥(¢) is continuous when t # 0. 
Then, for any set of real numbers, 21, 22, --- , 21, of which none is zero, 
there exists a positive number a, as scale, such that the likelihood condition 


(20) doi (a./a) ¢'(ai/a) + n = 0 
is satisfied. . 
The conclusion is also valid, if in place of the limit 8, there is postulated 


(27) lim y(t) = —a|o!]= lim y(d), 


t——b+0 t—>c—0 





where b > 0, ¢ > 0, and y¥(é) is continuous for —b < t < Oand for0 <t <e. 
That is, the new limits are to be infinite with sign opposite to that of a. 

Proof. ‘The limits for ¢ — 0 and for |t| — « are the same as the limits 
fora — « and a — 0+,—noting that t = x/a,x # 0. Thus 2y(t;) changes 
sign as a goes from 0+ to ©. Hence, since y(t) is continuous, (20) is satisfied 
for some positive a. 

For the proof of the second part of the theorem, suppose that x, > 0 and 
that z, is the greatest z;. Then witha > z,/c, but approaching x,/c, ¥(z,/a) 
becomes infinite with sign opposite to that of a. Furthermore, in 2y(x;/a), 
the positive x’s < zx, have a negligible effect; and thus lim y(z;/a), as 
a — (z,/c) + 0, is infinite with sign opposite to that of a, when this sum 2 
is taken for the positive z’s. Likewise, if 7; < 0, and is the least 2; , lim Ly(a;/a), 
as a — (—2,/b) + 0, is infinite with sign opposite to that of a, when this sum 
is taken for the negative x’s. If, now, the measurements happen to be all 
positive, we think of a as approaching z,/c + 0; and the continuity condition 
leads to an a which makes Zy(x;/a) = 0. Likewise, if the measurements 
happen to be all negative, we use —2,/b + 0. If both positive and negative 
x’s appear, we use the greater of the two ratios —2,/b and z,/c. 





















SECTION 3. SOME FAIRLY REGULAR FREQUENCY FUNCTIONS 





To illustrate the foregoing theorems in a somewhat general manner, consider 







the measurements, 21, 22, ---, %,, and with ¢ = z/a, t; = z,/a, set up the 
function: 
(28) f(t) = Ca”? | kt |? (1 4 kt?) e7 | ke . 


where, as before, C is a suitably chosen constant. 
Suppose also that 


(29) 









p> —l, q 20, 


and that either 
(30) 
Then with ¢(t) = log f(t), it follows that, when t ¥ 0, 

(31) to'(t) + 1 = (p+ 1) — rsk* |t |? — 2gk’(L + RE)”. 








r>G@er>s or r=0,2g>pt+l. 
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Now the condition 1-¢’(1) = —1 would be satisfied if ¥(k) = 0, where 
(32) W(k) = rsk’*? + rsk* + (2g — p — 1k — (p +1). 


But, under the conditions (29) and (30) ¥(0) < 0, and ¥(“) > 0. Hence, 
there is a positive k for which ¥(k) = 0. Then if k be assigned this value, 
(19) is satisfied; and by Theorem I, any scale a that the likelihood condition 
(20) may lead to isa mean. But, by Theorem III a scale a will actually exist 
—indeed, for any positive k that may be used in (29); since the limit of tg’(t) + 1 
is positive as ¢t — 0, and is negative as |t| — o. 

Moreover, if in (29), the further condition —1 < p S 0 is introduced, (22) is 
satisfied. And, thus,a < Maximum |z;|. Also, | é¢’(t) | increases with | t |. 
Hence, by (24) et seq., Minimum | z;| S a. 

If in (28), we set g = 0, s = 1, r > 0, and confine our attention to positive 
x and t, there is obtained the Pearson Type III. Reference to (32) shows that 
v(k) = Oif k = (p + 1)/r. With this substitution, 


(33) fhu  <*tcr™., C’ = Constant. 


Since ¢’(1) = —1, any solution of the likelihood condition is a mean. Here, 
with t > 0, t¢’(t) = p — (p + 1)t, and #¢’(t) — 1 = —(p +1). From (14) 
we see that, with p + 1 > 0, any mean obtained corresponds to maximum likeli- 
hood and the single maximum found is actually the largest value. Moreover, 
with the measurements, 27, 22, --- , Xn, all positive, a scale a will exist,—as 
noted in the general case (28). 

In passing, it may be noted that Type III appears’ rather naturally in a 


form giving ¢’(1) = —1 at once, without any transformation. Here, then, a 
scale is a mean. 


Given the Pearson Type I in the form 


(34) f(t) = Ca*(b + kt)"(c — kt)’, t=2/a, b>O0, c>O0, |pq|>O0. 


Ifp +q+1 > 0, it is possible to find a positive k so that with @ = log f, 
¢’(1) = —1. In this case, any scale found by the likelihood condition is a 
mean. With k thus chosen, f(t) has essentially the same form as it would have 
ifk = 1. Hence for convenience, let us simply set k = 1 in the above equation. 
Then for —b < t <¢, 


V(t) = to'(t) +1 = 1+ pttb+t)™ — ge — 1)”. 


Suppose now that p > Oandq > 0. Then Theorem III may be applied; since 
lim y(t) = 1, ast — 0; but lim Y(t) —- —~,ast— —b+0, orast—c — 0. 


6 Carver, H. C., Handbook of Mathematical Statistics, Chap. VII, see p. 105, Line 4, 
noting that ¢’ = y’/y. 
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Hence a scale a satisfying the likelihood condition exists. : Moreover, the likeli- 
hood is at a maximum; since, with —b < t <¢, 


o"(t) — 1 = —p'(b+ t)* — ge - )* - 1 <0. 


This maximum is also the largest value for all values of a. 
If the Pearson Type IV is given in the form 


(35) ft) = Cal + RP)? et ren, t= 2/a 


then if p > 1/2, it is possible to find a positive k which will make ¢’(1) = —1. 
In this case, any scale a is a mean. Moreover—for any k # 0—the limit of 
to’(t) + 1 is 1 fort > 0 andis 1 — 2pfort— «. Hence, by Theorem III, 
if p > 1/2, as above, then a scale a exists satisfying the likelihood condition (20). 


SECTION 4. FREQUENCY FUNCTIONS WITH CERTAIN PECULIARITIES 


The theorems of section 2 give sufficient conditions, which in some cases 
may not be necessary. Nevertheless, by violating certain hypotheses, particu- 
lar functions may be set up which exhibit various peculiarities. 

For the Pearson Types, the differential equation is 


, ao + at 
36 = GO « 4, t = 2/a. 
(36) vO bo + bi + det” 


/ 


The determination of a positive scale a by the Fisher likelihood process is 
impossible here, in case do = 0, a1 > 0, bo + bit + bet? > 0. For in this case 
to’(t) = 0; and thus (20) cannot be satisfied. The U-shaped Type II curves 
are in this class. Likewise, if a7) ~ 0, a, = 0, and bo + byt + bot” > 0,—for 
example, with bs > 0, bi < 4bob2 ,—and the measurements all happen to have 
the same sign as a,, such scaling is impossible. 

For the purpose of constructing peculiar functions we may take c > 0 and 
require that the measurements 2; be either —c or c—with at least one —c and at 
least one c—and that ¢(t) be an even function. Then ¢(—c) = ¢(c) and (11) 
becomes 


(37) Z on [Ca pr. 
The likelihood condition (13) reduces to 
(38) 0 = v(t) = ie’(t) +1 = (c/a)o’(c/a) +1, 


with the right member an even function of c/a. And from (14), a maximum 
likelihood is indicated when 


(39) (c/a)’$""(c/a) — 1 < 0, 


with the left member likewise an even function. A minimum likelihood is 
indicated if the left member is positive. 
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Let us apply this to the case where 
(40) $(t) = (—2/3) log (1 — 3/¢|);  #’@) = 2] t| (1 — 3] e])”. 


The likelihood condition (38) is satisfied only whent = +1. Also¢g’(1) = —1. 
Thus the only means are the internal means +c; and the only scale conformable 
to (38) isa =c. But this has minimum likelihood; since 1-¢’"(1) — 1 = 3 > 0. 
For positive ¢, this function (40) is a Pearson Type. 

Consider next a function of the form (28),—with p = —1.25, q = —0.5, 
however,—for which (31) becomes 


(41) @#(t)+1=-1/4-%/44+?/14+) = -a—- 9/40 4%). 


whence ¢’(1) = —1, ¢’(1) = +1, ¢’’(1) = —3. Here the likelihood condi- 
tion (38) has but a single absolute solution | | = 1, leading to the single scale 
a = c, and to the two internal means, +c. But, in this case 1-9’(1) — 1 = 0, 
so that a log L/aa” = 0. Moreover, fort = 1, a log L/aa®° = a * #0. Thus, 
the only scale obtained by the likelihood method (38)—viz., a = c— has a 
likelihood which is neither at a maximum nor at a minimum. 

Another anomalous function is that given by 


(42) o(t) = t* — 2.58, t = +c/a. 
The likelihood condition (38) leads to 
yt) = (1 —#)(1 — 4) = 0. 


The only solutions are t = +1, giving internal means +c; and t = + 1/2, giving 
external means +2c. And from (39) et seq., it can be shown that the internal 
mean and scale, a = c has minimum likelihood, while the external mean and 
scale, a = 2c, has maximum likelihood. 

But it will be noted that a maximum value for a vicinity does not always 
signify a largest value jor the entire possible range. Indeed, for the function 
(42), a = 2c has maximum likelihood without having the largest likelihood. 
To avoid such an anomaly, a necessary condition is that as |t| > ©, 
y(t) — — ©;as seen by taking the logarithm of L in (37), noting that asa — 0, 
(—log a) ~ +o. 

Finally avoiding the anomaly just mentioned, let us set up a frequency 
function, using the y(t) in (38), and writing 


y(t) =1+ (2) = (1 — 2°) — HU — 0.90). 
From this it follows readily that 


(43) o(t) = K — 1.95t? + 1.175t* — 0.3¢°, K = Constant. 


* 


This, with ¢; = -+-c/a, leads to an internal mean or scale a = c with minimum 
likelihood, a nearby scale a = c +/0.9 with maximum likelihood—differing 
indeed only slightly from the minimum just mentioned—and another scale 
a = cv/2 having maximum likelihood, and this likelihood is indeed greater 
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than that for any other positive value of a. The external mean a = cv/2 
in this case has the largest likelihood. This may be checked by the use of 
the logarithm of L as it appears in (37), in which the important part is 
(c/a) — log a. 

In passing it may be noted that if y(t) has the form y(t) = (1 — #)H(#), with 
H(1) # o, and t; = 2;/a; then any solution a of the likelihood condition 
y(t) = 0 is a mean,—by Theorem I. 


SECTION 5. SUMMARY 


When the R. A. Fisher likelihood method is used to find an “‘optimum”’ scale 
for frequency functions, it sometimes happens that this scale is a well known 
mean or at least is a substitutive mean—See Equation (5). Or a simple trans- 
formation (15) may often put the frequency function into such a form. Con- 
ditions are given under which a scale will be a mean. Under further condi- 
tions this mean will be internal—at least as regards absolute values. Finally, 
under certain conditions, a scale will exist. 

But for certain functions not satisfying these conditions, anomalies appear. 
The scale given by the usual likelihood condition may be a scale with a minimum 
likelihood. Sometimes the likelihood will be at neither a maximum nor a 
minimum. In certain simple cases, no scale exists. Furthermore, it may 
happen that the scales which are internal means have minimum likelihood and 
those that are external means have maximum likelihood. Among Pearson 
Types are found both anomalous functions and functions which would be 
regarded as regular as regards maximum likelihood. 

In this problem of scaling, likelihood is proportional to a posteriori probability 
with the a priori probability taken as constant. 








MOMENTS OF ANY RATIONAL INTEGRAL ISOBARIC SAMPLE 
MOMENT FUNCTION 


By Pau. 8S. Dwyer 


Introduction 


The problem of moments of moments has been investigated by a number of 
authors. The assumption of an infinite universe (or that of a finite universe 
with replacements) permits the application of the ‘‘algebraic’”’? method, the 
method of semi-invariants as introduced by Thiele (1) and developed by C. C. 
Craig (2) and the combinatorial analysis method introduced by R. A. Fisher (3) 
and used by N. St. Georgescu (4). A combinatorial analysis method has the 
particular advantage that it enables one to compute separate terms of a given 
formula. 

The formulae for moments of moments have been simplified through the 
use of new moment functions. Thiele introduced the half-invariant (1) which 
resulted in considerable condensation. More recently Prof. R. A. Fisher (3) 
has introduced the sample function k whose expected value is a half invariant. 
The most compact formulization presented thus far is his formulation of the 
half invariants of the sample k, in terms of the half invariants of the universe. 
This very compactness, however, makes it difficult to compare results with 
those expressed in the more conventional sample functions. Dr. Wishart has 
written a paper (7) in which he shows, among other things, how the Fisher results 
can be translated to the more conventional (Craig) results and vice versa, but 
such translation is in general no simple matter. It appears that the Fisher 
results are not immediately useful to the statistician who desires the formulae 
to be expressed in terms of the usual sample moment function. On the other 
hand the Fisher formulization is a remarkable discovery toward that harmony 
which must be naturally inherent in the field of moments of moments. Soper 
(6, 111) expressed the general situation when he wrote, “If the terrifying over- 
growth of algebraic formulation accompanying this branch of statistical inquiry 
is destined to have a chief utility in induction and going back to causes, then 
perhaps Dr. Fisher’s way of estimating a sample will prove to be most fertile, 
but if it is to be applied to problems of deduction, say to problems of suc- 
cessive eventuation such as propagation, then Mr. Craig’s plain moments seem 
to have a firmer hold on the exigencies of time.” 

It would appear then that the Fisher formulae and the Craig formulae are 
both needed. Georgescu (4) showed a partial connection between them in 
applying to the m functions a combinatory analysis somewhat similar to that 
applied by R. A. Fisher to the k function. It is the purpose of the present 
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paper to work out a combinatorial procedure for a more general sample function 
so that either the Fisher or Georgescu combinatorial results come out as special 
cases. In making such a generalization no limitation is placed on the sample 
function except that it be rational integral and that all terms are of the same 
weight. Thus the results are applicable to m,, m, + k,, m,k,, ete. as well 
mM, 
kk, 
the important formulae for the moments of a new sample moment function 
will be available by simple substitution as soon as any such new function is 
defined by a rational integral isobaric expansion of power sums. 

It is thus the purpose of this paper to determine the moments of a general 
moment function of the sample. This is done by keeping the multipliers of 
the various partitions of power sums indefinite until all manipulation is complete. 
It is then possible to assign the definite values of these multipliers which are 
associated with the desired sample function and to obtain the moment of 
the desired moment function in this way. Thus the Fisher result «(42), and 
the Craig result Si:(vs, v2) are special cases of the new result Au(fs, fe). It 
is obvious that it is not possible to carry the results using these general moment 
functions as far as Fisher and Wishart (3), (5), (7), have carried the results of 
the decidedly advantageous (from the standpoint of simplicity of result) k func- 
tion and yet it is surprising to find the simplicity which can be obtained in 
the general case. Incidentally the introduction of the more general symbols 
clarifies the successive steps of the partition analysis which are somewhat con- 
fusing in any specific case because of the insertion of the value of the coeffi- 
cients of the power sums in which the sample moment function is expressed. 

This paper is divided into three parts. The first part includes the necessary 
definitions, the basic formulae, and the general development of the algebraic 
method. In order to facilitate the algebraic work there is inserted a table giving 
the expected values of all possible partition products of power sums whose 
weight <8. The second part deals with the different sample functions which 
might be used. The third part gives a list of the various partition formulae, 
of weight <8, which contain no unit parts and shows how these can be used in 
writing the chief variations of the formulae for moments of moments. 


as to m, and k, although they are not applicable to ~/m, or 


In this way 


Part I 


1. General Moment Functions. Different moment functions have been de- 
fined in various ways, but all moment functions have in common the property 
that they may be expressed in terms of the power sums. It appears sensible 
to use this expression in terms of power sums as the working algebraic definition 
of moment functions. For example the function k;, which is defined by R. A. 
Fisher to be that function of the sample whose expected value is the third 
cumulant (half invariant) is to be given the working definition of 


m = — 8) _ QM), A) 
-" n—1DH—2)) (n— Dim — 2) n(n — In — 2) 
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where the numerical expressions in parentheses indicate power sums of the 
sample. 

“very term in the definition of a sample function has a “weight’’ which is 
equal to the sum of the power sums whose product is indicated by the term. 
Thus the weight of each of the terms of k; is 3. If all the terms of a given 
moment function have the same weight, the function is called isobaric and 
the weight of the function is equal to the weight of each term. Thus fk; is an 
isobaric moment function and its weight is 3. Since all the functions so far 
proposed are isobaric we limit this generalization of moment functions to iso- 
baric moment functions although it is possible that a more complex analysis 
could be worked out for non-isobaric functions. 

Generality demands the inclusion of every possible partition product of 
power sums. Such generality can be obtained by writing 


fi = a,(1) 
fo = a,(2) + ay(1)* 
fs = a3(3) + az(2)(1) + ain(1)° 
fa = a4(4) + a31(3)(1) + ag2(2)* + agi2(2)(1)* + ans(1)* 
and in general 
fe = Qo ants ok ots (pi)™ (pa)"* +++ (ps)"" 


where (p:)"' (pe)? --- (p,)"* indicates any partition product of power sums, 
@,%1 ... p%» is its coefficient and the summation is taken for every possible parti- 
tion. The number of parts of the partition is p = =r. It may be assumed, 
without loss of generality, that the partition is ordered, i.e. 


DP 2 pe = Ps 2--+ Z2Ds. 


A natural numerical coefficient of each term is the number of ways the r 
units can be collected to form the given partition. This value is given by 


r! 


1" 
® pi? --- i ~ (pi!)™ (pe!) «+» (pe) mi! me! +» 


If we set 
Tr 
ae T l TT Tr 
Api)... ps = ™ ee Apil.-- pet 
Pi ece P. 


the definition of f, becomes 


r= D( 


In the present paper the capital letters are used to represent the corresponding 





* 


) ay. pte (pi)"* +++ (ps)"* 


T T 
pi’ ++: Ds! 
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functions of the universe as defined by the corresponding power sums of the 
universe. Thus 


i‘ 


Ww 


Pi -2+ Dp," 


) ayn ge GP" --- Fy” 


represents the corresponding function of the universe. In the case of the 
moment about the mean and the semi-invariant the Greek letters » and \ have 
been used to represent the corresponding function of the universe. In the 
case of functions whose notation is quite widely established, it is preferable to 
use the conventional notation, but in introducing new functions it appears 
wise to use the relationship between small and capital letters since the corre- 
spondence between the English and Greek alphabets is not exactly one to one. 
It should be particularly noticed that this notation does not agree with a pre- 
viously accepted scheme of using the small English letter to indicate the function 
whose expected value is indicated by the corresponding Greek letter. In the 
present paper it is not the expected value property which serves as the basis 
of notation but rather the definition of the function in terms of the partition 
products of power sums. 


2. The Working Definition of Moments About a Fixed Point. The sample 
functions defined by 


» (1) » (2) , (3) » (r) 
m, = —, mM, = —, My = —,-+-, Mm, = — 
n n n n 


are obtained from f, by placing 


(1 
— when s = 1, m = 1, and p; =r. 
in 


Ls =—_ 
-- pst = 


0 in all other cases. 
The Greek y’ is used to indicate the corresponding function of the universe. 


3. The Working Definition of Moments About the Mean. The moments 
about the mean are defined by 


r (1) 
m = —, Mme = — ——) 


n n 


Ms 


9 


n n? n>’ 


_ 3) _ 32)0) , 2a" _ ) _ 4) 
n ? 


n- 
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and in general m, is obtained from f, by placing 
(1 
in 
(— tan >1,m = 1,8 = 2, and pz = 1. 


1+ 
én... <\|\om 


i¢-iF" © -® 


- —ifp, = 1,s = 1,and™ =r. 
n 


ifs = 1,7, = 1, andp: =r. 


(0 in all other cases. 


The corresponding moments of the universe are indicated by the conventional yz. 
For conciseness moments about the mean are referred to as ‘“‘moments.” 


4. The Working Definition of the Half Invariants. The half invariant 
moment functions of Thiele, as applied to the sample power sums are [see C. C. 
Craig (2, 7-10) and Frisch (12, 20-21)] 


yD @) _ 


L = — a , = &) — 32) 0) , 20) 


9 


n’ n n? n n? n3 


1, — 4) — 4) G@) _ 32)" | 12(2) 4)" __ 6G)" 
pe st oS SS 


n n? n> ni n4 


and in general 


S| I" Siac . 
l, = (pi)"* (pe)"* «++ (ps) 
Zi n? pPi' +++ pe’ Ps) P P 
so that 
_ (-1""'@- 0! 

n? 


a 


pil --- p38 
The corresponding moments of the universe are indicated, after Thiele (1) 
and Craig (2), by \. R. A. Fisher (3) used x while Georgescu (4) used s. 

In the present paper these functions are referred to as ‘Thiele moments.” 


5. The k Functions of R. A. Fisher. The k statistics of R. A. Fisher are 
defined in terms of the sample power sums by 
ee 
aa .  f-—) o-oo 
», — 23) _ 82) ) ss, 2)" 
*“(n-DY(n—-2) (n—1(n—-2)° n® 


k, = MVE) _ 4m + YE) 3(2)° 12(2) (1)* __ 6(1)' 


(n — 1)® m—-D® KM —2°9*M—1)% no" 
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These values and values for ks and ke are given by R. A. Fisher (8, 203-4) 
while algebraic methods of attaining them are presented in sections 16, 17. 
They are referred to as Fisher moments. The corresponding functions of the 
universe, if used, would be represented by K,. 


6. The A Function. Just as Fisher introduced a sample function whose 
expected value is a Thiele moment of the universe, so it is possible to introduce 
a function whose expected value is a moment of the universe. Such a function 
is defined by 


,_ (1) (2) (a) 
on io <i 
h n(3) 3(2) (1) ‘ 2(1)° 


sina (n — 1) (n — 2) os (n — 1)(n — 2) n®) 
he = (n® — 2n + 3) (4) _ 4(n° — 2n + 3) (3) (1) _ 3(2n — 3) (2) 








(n 1)® ; n4) ; n 
6(2) (1)° 3(1)* 
(n—-1)® 7" 


Methods of obtaining the expansion of this function in terms of power sums 
are presented in section 18. The corresponding function of the universe, if it 
were used, would be represented by H,. 


7. Other Moment Functions. It is possible to obtain an indefinite number of 
moment functions. For example one might define a function of weight 2 whose 
variance equals yy, (or u2). It is possible by the methods of this paper to 
find expressions for such moments. 

For reference purposes Table I is provided showing the values of a for each 
partition of weight <6 for the functions m’, m, l, h, k. The values of 


( : 
pi’ po’ +--+ p,* 


are also inserted, in the left hand column, so that it is possible to read from the 
table the values for f = m,, m,, l,, ky when r < 6. 


8. Products of f Functions. The product of two or more isobaric functions 


is also isobaric and of weight equal to the sum of the weights of the functions. 
Thus 


fof = [a2(2) + an(1)(1)J[a(1)] = a2ai(2)(1) + anai(1)° 
fofi = aya;(2)(1)° + aya\(1)*. 


In multiplying f,, by f,, any term of f,, is of weight 7, and when it is multi- 
plied by any term of weight 72, the result is a term of weight 7 + re. 
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TABLE I 
Coefficients of Products of Power Sums in the Expansion of Different Moment 


Functions 


k, 
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n(n + 1) 
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TABLE I—Concluded 


| 
Numeri- 
cal 
coeffi- 
cient 





2(n + 2) 
” 7 (1 


2(n — 1) 
in — 
— 6 —_— 
24 4 


n) | n®) 





15 








2111 


Q@j1111 











R. A. Fisher [3, 207] used the product kk: as an illustration of the algebraic 
method. The more general f3fo gives 


Safe = [a3(3) + 3a(2)(1) + ay(1)*}'[a2(2) + an(1)(Q1)] 

= a5ax(3)(3)(2) + aban(3)(3)(1)(1) + Gasaaar(3)(2)°(1) 

+ [6asa2ay, + 2asded1](3)(2)(1)* + 9a3,a2(2)°(1)* + 2asaya11(3)(1)” 

+ [6a2d:n42 + 9a21an](2)°(1)* + [6a2ainan + a2ain)(2)(1)® + ainan(1)* 
which reduces to the value as given by him when the values of a are substituted 
from Table I. 

9. The Expected Value of Any Partition Product. The expected values of 
partition products are well known and are indicated by 
E(~1) = Nip 
E(p:)(p2) = Mtp.+p2 + (mn — Wpp Hp» 
E(p:)(p2)(ps) = NU py+02403 + (NM — 1) [up 42H ps + Mortrattns + Moot rattnrl 
+ n(n — 1) (n — 2) wh wpottys- 
and in general 


r eee Ds 3 
E(p:)"*(p2)"* «++ (ps)"* = Zz n?/ 2 xe 


-+ Gi 


Ts 


(ul, Gl oo Gul 


Ts 


8 
indicates the 
x X2 x 
qt’ qz* +++ git 


pi' po? ---p 


where +t = x1 + X2 + X3 +-°- + Xe and ( 
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number of ways in which the partition pi' pz? --- py* can be grouped to 
form the partition gf' qx? --- qx'. 

The continued application of the result above leads to a large number of 
formulae. In order to make these results accessible I present in Table II the 
expected values of all partition products of weight <8. The essence of the 

T1 2 
table is the evaluation of the expression _ . The numbers 
qi' qx? --- qi' 
at the top of each column indicate the subscripts of the u’s which must, of 
course, be multiplied by n. The entries on the extreme left are the numerical 
coefficients associated with each row. 


10. The Expected Values of the f Functions. With the use of Table II one 
is able to write expressions for the expected values of f, when r < 9. 


w(t) = E(fi) = ana; 

ui(fe) = E(fo) = (a2 + an)nys + aun(n — 1)uy? 

ui(fs) = E(fs) = (a3 + 3ax + ain)nys + 3(ae1 + ain)n(n — 1)pomi 
+ ayn(n — 1)(n — 2)y;’ ete. 


If the expected values of the f functions are expressed in terms of the moments 
about the mean of the universe, these formulae become, since u; = 0 


ui(fi) = 0 

ui(fo) = (a2 + au)rms 

ui(fs) = (as + 3an + ain)nus 

milfs) = (as + 4a + 3a22 + 6a2n + aun)nus 


+ 3(a2 + 2a + aun)n(n — 1)ys ete. 
These may be written more symbolically as 
ui(fi) = 0 
ui(fe) = bonus 
ui(fs) = bars 
ui(fs) = byrpia + Bdoon(n — 1)ys ete. 


11. The Expected Value of Products of f Functions. The expected value of 
products of f functions may be similarly found. For example 


u2(fe) = E(f2) = Elae(2) + an(1)*) = a2H(2)” + 2aanE(2)(1)(1) + aiE(1)*. 










TABLE II 
Expected Values of Partition Products 








z ~- 
coef. | nm | n&2)) (3)| 
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x | 3/21/2111 


























1 | 3 
|} 3 | 21 
1 | 2B 
weight = 4 
jeoef.| | m | n® n®| n®) ne 
| x | 4 | 31 | 22/211) 18 | 





























weight = 5 
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: ~ 
| coef. | n | n®) n®) n@) n® | nM) n(5) 









7 | 5 | 41 
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32 | 31? | 221 | 213] 15 | 


























weight = 6 
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TABLE II—Continued 






















































| coef. | n | n®) n@) n®| ni n®| ns) ni*| n(4) | ni5)| ni6) 
| | | 6 | 51 | 42 | 33 411/321) 2° | 313| 22 12] 214) 18 
6 | 51 

15 | 42) 

10 | 33 

15 | 411 

60 | 321 

15 | 23 

20 | 31° 

45 | 2? 12 

15 | 214 

1 | 18 


weight = 7 






























































































| coef. n | n)| ni?) nia) ni) ni®| ni) | n(3)| n| ns) | m4) | 6) | n(5) { nie ni| 
| | 17 61 | 52 | 43 | 512 | 421| 32 1 | 322 419 | 3212 | 29 1 | 314 | 22 13 21s | 17 | 
| 1 |7 | | 
7 «| 61 
21 | 52 
| 35 | 43 
/ 21 | 512 
105 | 421 
| 70 | 321 
| 105 | 32? 
135 | 41 
210 | 3212 
105 | 231, 
135 lave fal 4 
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Table II can now be used by indicating a3 as a multiplier of E(2)’, 2a,ay, as a 
multiplier of E(2)(1)(1) and aj; as a multiplier of (1)*. Then at once it is 
evident that 

u2( fo) = (a) + 2acay, + ay) Nu4 + (a3 + 2acay, + 3aj1)n(n — 1) us 
(a2 + an)°*nus + [(a2 + an)” + 2an)n(n — 1)y5 
bonus + (b2 + 2bii)n(m — 1)ys. 


Similarly 
un(fs, fo) = bsbonus + (b3b2 + 3beibe + 6babu)n(n — 1)psue 


us(fs) = bsnus + (9ba1 + Gbaber)n(n — 1)usus + (bs + 9bi)n(n — 1)u5 
+ (9b3, + 6bin)n(n — 1)(n — 2)u3 


etc. 
where bs = 3 + 3d2 — ai, bey = Ag) 4. ai, bin = A111. The important 


special cases are obtained by assigning the proper values to the a’s as given 


in Table I. Thus 
’ 1 
u2(M2) = 73 [(m — 1)?us + (n® — 2n + 3) (n — 1)y5] 


which agrees with the corrected result of “Student” in 1908 (8, 3) and Tchou- 
proff (10, 192). Similarly 


ui1(M3, M2) = Zt — 1)*(n — 2)us + (n — 1) (n — 2) (n? — 5m + 10) psue] 


us(ms) = 5 [(n — 1)°(m — 2)°us + (—6n + 15) (n — 1) (n — 2)° sue 


+ (n® — 2n + 10) (nm — 1) (n — 2)?n5 + (On? — 36n + 60) (n — 1) (n — 2)yi] 
ete. 

In the same way 

us , (n> — 2n + 3)mr 


ua(ke) n n(n — 1) 














‘77, " Ms (n? — 5n + 10) uspe 
wiulks, ke) 7. n(n — 1) 


us(ks) _ 6 4 (—6n + 15)psue (n’ —2n + 10) u3 (On? — 36n + 60)u3 
7 








u n(n — 1) + n(n — 1) “n(n — 1) (n — 2) 


ete. 
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and 


us(ms) = ~ we + (x — Dud 
, , / 1 
Mii(m3, Me) = - [us + (nm — 1) papel 
, , 1 : 
uo(m3) = ra [us + (n — 1)ys] 


ete. 


12. The Expected Value of the Products of f Functions in Terms of the 
Thiele Moments of the Universe. The formulae giving the y’s in term if the 


rm gg 


d’s are 
Me = de 
M3 = Xz 
va = Ag + 33 


Ms = As + 10A3A2 
us = Ae + 15Adde + 103 + 15d2 


1 
— >( ai .) (dpi)? py)" «++ Op.) 
Pi P2 o&'s Ds 


where the summation holds for those partitions having no unit parts. See 
the results of Craig (2, 7-11) and Frisch (12, 21). It is at once possible to 
express the moment formulae in terms of the Thiele moments of the universe. 
Thus the general results above become 


u2(fo) = bands + [Bban + (b2 + 2bi)n(m — 1)]r2 
un(fs, fo) = dsbends + [10bsbon + (babe + 3baib2 + 6babu)n(n — 1)]Asdz 
us(fs) = b3nde + [15b3n + (9bir + 6bsbn)n(n — 1)]Aadre 
+ [10bsn + (63 + 9ba)n(n — 1)]d3 
+ [15b3n + (27b3; + 18bsb21)n(n — 1) + (90a: + 6bin)n(n — 1)(n — 2)JA3. 
13. The Thiele Moments of the f’s in terms of Thiele Moments. It is 


now possible to reduce to the Thiele moments of the f’s by means of the usual 
relations 


do(f,) = ua(fr) — wi (S-) 
An(fr; ’ Sra) = wifey ’ i) ram mio Sr, ’ Seo) mor Sr ’ >) 
Aa(f-) = us(f-) — 3u2(f)ui(f-) + 2mr°(f) 
etc. 





” EE gg 
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so that the results become 
Ao(fo) = bands + 2[ban + birn(n — 1)]r2 
Au(fs, fo) = bsbondrs + {3[bsbon + boiben(n — 1)] + 6[bsbon + barbun(n —1)]}rsrz 
No(fa) = bands + {6[b3n + bsbarn(n — 1)] + Y[bsn + bain(n — 1)]}rsdrz 
+ Q[b5n + biin(n — 1))A3 + {9[bin + Qsban(n — 1) + biin(n — 1) + ban™] 
+ 6[b3n + 3ban(n — 1) + biun(n — 1)(n — 2)]}r2 
etc. 


The formulae as written are adapted to the partition representation of Part ITI. 
When the f’s are equal to the m’s we have 
(n—1) , An— Drm 
n> v n2 - 


(n — 1)?(n — 2)As + 6(n — 1) (n — 2)Asre 


nt n 


Ao(me) = 





Au(mMs, mM) = 
2 9 2 
dems) = (n — 1) : — 2). + 9(n — = 2)*rsre 
n n 
+ 9(n — 1) . — 2)°r5 ‘ 6(n — 1) (n — 2)r3 
n n 
ete. 


which are the results as previously given by C. C. Craig (2,55). In like manner 
when the f, = k, 











Na 2r3 
Nelle) = = + ——— 
rn(k we 6 AsA2 
ew  . a | 
— ir 9x3 6n A> 
Whe) = + its -ita-pe-® 


ete. 


as given by R. A. Fisher [3, 210] while 


, 1 
Ao(ms) = 7M + 2n3) 
, , 1 
Au(m3, M2) = 7 (ns + Qsdz) 


peleat) an - (ng + 15dsde + 9A2 + 15r3). 


ete. 
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14. Various Formulization of Results. Although different moment functions 
of the universe may be used it is customary to express the results in terms of 
universe moments about a fixed point, in terms of universe moments, or in 
terms of universe Thiele moments. It is possible to express results in any of 
the nine forms 


u’(f,) (moments about a fixed point (u’) 


u(f-) > in terms of ae (u) 


(fr) ) 
where f, represents the isobaric sample moment function of weight r. One 
purp%se of such varied formulization is to discover the most compact form 
and also the one best adapted to use in the case of a normal universe or a uni- 
verse whose moments obey some discoverable law. As suggested above Craig 
(2) has shown the relative compactness obtained by using A(m,) and Thiele 
moments of the universe while R. A. Fisher (3) has shown the great additional 
compactness obtained by taking f, = k,. 


Thiele moments (A) 


15. The Application of the Algebraic Method to ai(f;, fe). Before leaving 
the algebraic method it is perhaps wise to outline the steps in the case of a 
more involved problem. We take the example which R. A. Fisher (3, 207) 
has used in the case in which f, = k,. To find Ag(fs, fe). 

The value of f3f2 was found in section 8. To find its expected value it is 
only necessary to enter the coefficients of the different partition products in 
this expansion at the left of the corresponding rows as indicated in Table II. 

The coefficient of any moment partition of the universe is found by multi- 
plying each column entry by its corresponding left row entry and then by 
multiplying by.n as indicated at the top. Thus the coefficient of ys is 


(a3az + aay + 6a3d2102 + 643021011 + 203011102 + Jasid2 + 2a30;1011 + 6a2dqdi1 
+ 9ae1010;; + 6421041101. + i 1102 + ajay) Nn 
which after some algebraic work reduces to 
(a3 + 3a + ain)’ (a2 + an)n = b3ben. 


In this manner it is possible to write the result either in terms of universe 
moments about a fixed point or in terms of universe moments. If moments 
are used, one may neglect all column partitions involving unity. 

It should be noted that the a’s defining k, as given in Table I can be inserted 
here if desired. If these multipliers are introduced throughout the rows and 
columnar partitions involving unit parts are not used one will arrive at Table I 
of R. A. Fisher [3, 208] though there are some slight typogranhical errors in 
his rows for (3)° (1)° and (3) (2°) (1). 

Determining all the coefficients in this manner we find after considerable 
algebraic manipulation that 
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ur(fs, fo) = bsbonus + [b3bo + Qboib2 + 12bebeibs + 6bsbarbeln(n — 1) sue 
+ [2b3b2 + 18baibe + 18boi:bi + 6bsbeibs + 12bsbabuln(n — 1)usus 
+ [2b3bu + Qboibe + 18baibu + Gbsbabeln(n — 1)us + [36b2ib2 
+ 54dsibi + 6bsbaibs + 12bsbabu + 12bsbinbu + 72barbinbu 
+ 18binbeln(n — 1)(n — 2)ugus + [D3be + 6bsboibe + 12bsbarbn 
+ 27berb2 + 9Obabu + 36barbinbe + 72babsubu + 36binbuln(n — 1)(n — 2)y5u2 





2 


+ [9boibe + 18b3:b11 + 36babiibu + 6bi11b2 + 36D} bi|n(n — 1)(n = 2)(n — 3) ua. 








If f, = k, the proper values of b are inserted and the expression above becomes 
that given by R. A. Fisher (3, 208). For example the coefficient of 2 is 











(Qn? — 63n” + 240n — 420) (n — 3) 
n(n — 1)?(n — 2) 





when 


1 1 1 
=> 6] ae es 














1 2 


ail n(n — 1)’ —- n(n — 1) (n — 2)’ 


The algebraic results involved in changing the general formula above to 
other functions are too extended to present here. A symbolic means of attaining 
them is included in later sections of the paper. 


Part II. The Determination of Specific f Functions 


16. Functions Determined by the b’s. In Part I it was shown how various f 
functions are defined by giving definite values to the coefficients of the power 
sums. It is the purpose of this part of the paper to show how functions can 
be specified by means of their expected values in terms of moments of the 
universe. This is essentially the method used by R. A. Fisher in defining his 
k function and it is here extended to other functions. In this case the b’s are 
first determined and the a’s are then found from them. The first moments 
of fi, fe, fs were given in section 10. To these we add, as shown by Table II 


ui(fs) = (as + 443; + 320 + 6aen + Ay) ny + 4(ds: + 3aen + aun)n(n — 1)usui 
+ 3(@2 + 2ae1 + aun)n(n — 1) us” + 6(aen + aiun)n(n — 1)(n — 2) uous” 
+ auun(n — 1)(n — 2)(n — 3) ui 


etc. 
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These can be written more symbolically in terms of the 6’s 

wi(fr) = binws 

ui(fe) = bonus + bun(n — 1)py 

ui(fs) = bsnus + 3ban(n — 1)wou, + birn(n — 1)(n — 2)u; 

wi(fs) = banws + 4bsin(n — 1)ugur + Sbeen(n — 1)us” + 6beun®? uour” + 6n uy*, 


and in general 


a 
, 7 t ° ’ 
u(y) — a = * bpm eee ps8 n(n.) (Ms)? ili iar”. 
Pi po +++ Ps’ 


The expansion of the function in terms of the power sums of the sample demands 
the determination of the a’s. This can be accomplished by solving the equations 


aq, = by 
dz + an = be 
an = bu 


a3 + 3d2 + Qin bs 


Ax + Qin = be 


Qin = bin 
dy + 443, + 3ae2 + 6en + Gun = da 
31 + 3de1 + din = bs 


dog + 2den + Gin = bee 


ete. 
The solutions are 
ay= b; 
dz = be — bu 
ay = bu 


ay = by — Shun + Som 

AQ = ba — bin 

Qu = buy 

ay = by — 4bn — Shee + 126 — Chin 
= bn — Sb + 2dnn 


bee = 2be11 + bun 


8 
| 


Q 
2 
3 


dor = ber — bun 


Qun = bun. 
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The values of a,, at least for r_< 4, follow the law 


«= D( 


r 


- . (—1)?" (p —1)! byt... pt 
P, ae D.° 


and 

da, = laza;! where \a,a;' indicates that a2 = bs — by is multiplied by a, = bi, 
the rule of multiplication being suffixing of subscripts. Similarly a2. = ‘asa! = 
(be — by) (b2 — bu)! = be — 2ben + bun. 

This statement illustrates a general theorem which will be established later 
in another paper by a different approach that for all cases. 


| 





r 


. . ) (—1)?" (p — 1)! bpm... pt 
pi’ oe ps’ 


and that 


Gr, ei ‘Gy, Gr, — a, , 
This theorem enables one to write, with comparative ease, the coefficient of 
any product of power sums in a sample function whose expected valued is defined. 


For example the functional coefficient of (3)(2) in fs is 





'azd2! = |(b3 — 3ba + 2bin) (bo — bu)! = 032 — bs — 3be01 + 5b — 2b 


while that of (3)(1)(1) is ‘agaya;) = bau. — 3be1n1 + 2binn. If the expected value 
of the function is known the b’s are determined and the values of the above 
expressions can be found by substitution. 


17. The Values of the Fisher Moments (k functions). The k functions have 


been defined to be these functions whose expected values are the Thiele moments 
' a , ‘ 
of the universe. Thus u:(k,) = A, and since 


= D( 


« . . / . . 
it follows at once that by comparison with yu; (f,) in the last section, that 


" 


rT 3 


T1 To 
Pi Po +--+ Ds 


)(-pe (o — 1)! (wy,)(upe)”? - +> (up,)™* 


(<=! 


bp pt? --. ps = — 


ni?) 





| 
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The insertion of these values in the formulae of section 16 gives the values of 
a such as those indicated in Table I and in section 5. Thus the coefficient of 


(3)(2) in fs is 
1 2 6 30 48 
10(b3e — bs ie 3be01 + 5bein ai Qbii1) ar 10 E + n®) n® 7 n) + 3 | 


— _ 10n® 
= (n _— 1)" 


The coefficient of (3)(1)(1) is 


48 10(2 4 
10(bsu — 3bon + 2bi11) = 10 E + = +> SI — = +2 


n n (n —_ 1) J 


18. The h Functions. It is also possible to define a function whose expected 
. . , 
value is the moment of the universe. Thus y:(h,) = yu, where 


“=D. 


pi’ ein Pp." 


/ , / 8 
Aggtr ... pte(ui,) (ng)? «+ * (up,)” 


(1ifs=1,™ = 1, and p; = r. 


| (—1)" if pr >1,m = 1,8 = 2and pp = 


wT 
o- ps* = 


(—1)"" (r — 1) if p, = 1,8 = 1, and m =r. 
\0 in all other cases. 
Comparing with the value of y;(f,) in section 16 we have 


A <r a 
by7 ings = aol welt 
nr) 


s 


The substitution of these values of b in the results of section 16 gives the expan- 
sions of h, in terms of power sums as illustrated by the formulae of section 6 
and Table I. Thus the coefficient of (3)(2) is 


10(b32 — bsi1 — 3be01 + Sbern1 — Zr) 


1 5 8 —10(n — 2) 
--wlo+ 5404545] - er? 


Similarly the coefficient of (3)(1)(1) in hs is 


10(n? — 
10(b3n — 3ben + 261m) = 10] ; + z + : | = = ae 2. 


n'3) n) n® mn 


19. The h’ Functions. One line of attack calls for the introduction of new 
moment functions which will result in simpler formulae. Thus for example, 
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C. C. Craig wrote (2, 37) “It rather seems that the best hopes of effectively 
further simplifying the problem of sampling for statistical characteristics lie 
either in the discovery of a new kind of symmetric function of all the observa- 
tions which may be used to characterize frequency functions and which will 
be more amenable than either moments or semi-invariants for use in sampling 
problems, or in, what may very well prove to be much better and more 
feasible, the abandonment of the method of characterizing frequency functions 
by symmetric functions of all the observations altogether.” 

R. A. Fisher has shown that it is possible to introduce symmetric functions 
which do simplify the resulting formula appreciably. It is the purpose of this 
section to introduce an additional symmetric function which simplifies the 
resulting formulae to a much greater extent. It is admitted that this function 
does not have all the properties (such as invariance with respect to change of 
origin) possessed by the Thiele and Fisher functions, but it does not have the 
property of making the resulting formulae simple. It also has the advantage 
that u(h,) = u’(h,). 

The basic idea is to find a sample moment function whose expected value is 0. 
A first attempt, placing every b = 0, is of no avail since every a is also equal 
to 0 and there is no function. A second attempt is based on the idea of finding 
the function h whose expected value is u;. If the universe is assumed to be 
measured about its mean, as is conventional, it follows at once that yu; = 0 
and ui(h,) = 0 so that 


ttge(Iie, » thes) = pye(Ir,» rs) 


This function then has the property that its moments about a fixed point and 
its moments are identical. 


In order to discover its expansion in terms of power sums, we note 
oe ty 
ui(h,) — on 


1 


n*) 


and 6,71... ,%s = 0 in all other cases. The a’s are determined in the usual 
way. Thus 


and it follows at once by comparison with y;(f,;) in section 16 that bi 


1 


ae n(n — 1) 


- 1 
~ n(n — 1) 


ll 


so that 


, 1 
ho = n(n — 1) [(2) - (1)(1)]. 
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Similarly 


, 


hs = — (2(3) — 3@2)(1) +" 


/ 


hy = as (6(4) — 8(3)(1) — 3(2)(2) + 6(2)(1)(1) — (1) 


and in general 


C.i2 1c 1 (mr — DN (ae — + Ke — DT 


n” 
1’ 
( Tr} .) (p.)" Eos -? 
Pi om D.* 


In order to show the simple form in which results can be given we substitute 
. . y / / 
the values of the 6’s in the results obtained above. Not only does y;(h,) = 0, 
but by section 11 


lt) = ph) = of = ——*—. of 
n(n — 1) 
An(hs, hy) = uulhs, hy) = wilhs, hy) = 
6 3 


do(h3) = us(hs) — uo(hs) = n(n — 1)(n — 2) Me 


while from section 15 


36 uzue 36(n — 3) us 











2 : hs = 2 : : = : 3 3 $$$ ——_—___—. 
Aai(hs, he) = war(hs, he) = waiths, he) n?(n — 1)? (n — 2) n(n — - 1) (n — 2) 


It is to be noticed that these formulae contain very few terms and that the 
terms themselves involve very low moments of the universe. This simplicity 
has been attained without making any assumption such as normality, regarding 
the nature of the universe. 


20. Table of Values of b for Different Functions When r < 6. This process 
of defining functions by means of expected values could be extended indefinitely. 
Perhaps it has been applied to enough functions to suggest the breadth of the 
applicability of the theory developed in Part I and Part III. 

As the b’s are the quantities which are used in the formulae I have provided 
Table III giving their values for the six functions, m}, m,, l,, k-, hr, he when 
r = 1, 2, 3, 4, 5. When the a’s are known, the b’s are computed from them 
according to the formulae of section 16. 
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TABLE III 
Values of the b’s forr < 5 
Num. b m f Mr ly ky hr h 
coef. r 
nate L 1 Lf} i fa 
n n nr n n n 
1 be 1 “2 o=% 1 1 0 
n n? n? n n 
“ta Pa 1 1 1}; 2/1 
t) ou | ° nt — ~ 2®| ~ n@| ne 
1 (n — 1) (n — 2) (n — 1) (n — 2) 1 1 
1| b& | = woe a ~ - | 0 
n n n n n 
(n — 2) n—2 1 1 
3| bu | 0 = ies i —7@|—pw| ° 
| 2 2 2 2 |1— 
1 bin 0 n3 n3 n® n® n® 
1 | (n—1) (nt —3n +3) | (n— 1) (n? — 6n +6) 1 1 
1 be a ee a —— SS a ai ms 0 
n n n n 
| (n? — 3n + 3) (n? — 6n + 6) 1 1 
4 ba 0 7 ns = ns Par ~~ may n®) 0 
2n -— 3 (n? — 4n + 6) 1 
3| be | 0 —— - =i @¢ | 
, n—3 2(n — 3) 2 1 
6 | bn | 0 = “a n® | n@ | 9 
3 6 6 1 
1 bin 0 ~ a! = ns n® | n™® n® 
vl) ge | 2 [= Wn = 2) (22m + 2)|(n—1) (n—2) (n?—12n-+12)} 1 .t. 
"is n> n> n n 
7 (n? — 4n? + 6n — 4) (n? — 14n? + 36n — 24) 1 1 
eee Se —— — n® | —~ p@ 0 
_ = 3 — Sn2 is 1 
10 | be | 0 e—2s2 a th oul © 0 
n n n 
: n? — 3n+4 2n? — 18n + 24 2 1 
10 bs 0 — a a n® n® 0 
| a J as 
| ats —_— 
| om lo,  —%>2 mime | 2) 4 | 
| n n n* 
n—4 6(n — 4) 6 1 
| | 4 84 24 4 1 
. oon) S ni ni n® | n® | nO 
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Part III. Combinatory Methods 


21. Partition Representation of Expected Value of f Functions. The formulae 
wi(fr) =: dinui 
ui(fo) = bonus + bun(n — 1)mr” 
ui(fs) = bsnus + 3ban(n — 1)usur + biun(n — 1)(n — 2)y;° 
ui(fs) = banws + 4ban(n — 1)uous + 3beaen(n — 1)u2” 
+ 6ben(n — 1)(n — 2) uous? > Dacsstt™” me" 
are “synthetically” given by the column partitions 


1 


2 


The partition parts represent both the subscripts of the moments and the 
subscripts of the b’s. If p indicates the number of parts, the n multiplier 
isn”. The numerical coefficient is obtained by taking the sum of the entries 
in the column (the weight) and dividing it by the factorials of all entries times 
the factorials of all repeated entries as indicated by 





1’ r! 
pips? + pre) (pi!) (pe!) +++ (ps!)"* mi! me! +++ re! 


The translation from the synthétic partition form to the expanded form is 
accelerated if the coefficients are known. These are provided in the following 
partition representation of the formula for 4;(f,) when r < 8 and the results 
are expressed in terms of the moments of the universe 


wa(fi): 
u3(fo): 


ui(fs): 
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wi(fi): 


ur(fs): 


ui(fe) : 


ui(fz): 


ui(fs): 


The proper formula can be stated immediately from its synthetic representa- 
tion. Thus for example 


ui (fe) = benus + 15byn(n — 1)ysue + 10b33n(n —.1)u3 
+ 15bo00n(n — 1)(n — 2)ud. 


22. Partition Representation of the Expected Value of a Product of f Func- 
tions. Two column partitions may be used similarly to represent the expected 
values of the products of two f’s, three column partitions for the expected value 
of the triple product, ete. In order to obtain all terms it is only necessary to 
combine every partition of each f in every possible way. The synthetic repre- 
sentation of E(mez, mj) is 


1 1 2 1 
21 20 11 10 
01 10 10 

01 


The sum of the entries in each row indicates the proper moment while the | 
number of rows indicates the number of parts as in the preceding section. 
The n coefficient associated with a p rowed partition is then n®. The b coeffi- 
cient is indicated by the columnar entries. Thus 


un(fe, fr = bebinust [bob, + 2bubiln(n — L)womr + bubm(n — 1)(n — 2)us*. 
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We verify this by the algebraic method 
uin(fe, fr) = Bf [ax(2) + an(1)(1)Ila(1)]} 
= E [aza;(2)(1) + anai(1)") 
= aea;[nus + n(n — 1)uoms 
+ ana[nus + 3n(n — 1)ysur + n(n — 1)(n — 2)y;°) 
(a2 + an)arnus + (a2 + an)an(n — 1)yom1 
+ anamou: + anan(n — 1)(n — 2)y;° 
= bebinus + bebin(n — 1)yaui + Qbubyn(m — 1)yomi 


+ bubm(n — 1)(n — 2)p;° 
as indicated. 


It thus appears that the partition representation is a mnemonic device for 
indicating the solution as obtained by the algebraic method. A more formal 
justification is based upon the property that if 


E(fe) = be(2) oe bu(1)(1) and E(fi) = b,(1) 


then E(fe, fi) can be obtained by a symbolic multiplication of be(2) + bu(1)(1) 
by 6,(1) where the b’s are multiplied but the power sums are collected in all 
possible ways. Thus 


E(fo, fi) = bebsl(3) + (2)(1)] + bubs[2(2)(1) + (1)] 


which gives 
E(fe, fi) — bebynus + bebin(n — 1) wou + 2biubin(n —_ 1) uous = bubsn®? us? 


as before. 

This symbolic multiplication is generally true and serves as the real algebraic 
justification of the partition representation. It will be established in a later 
paper dealing with the more general case of a finite population. The general 
type of partition analysis has been used previously by Fisher (3) and Georgescu 
(4). Each has established it through analytic rather than algebraic means. 


23. Determination of the Coefficients. Methods of determining the numerical 
coefficient have previously been given by such authors as Fisher (3), Wishart (5) 
(7) and Georgescu (4). If the f’s are of different weight, the coefficients of any 
partition (an interchange of rows is not looked upon as changing the partition) 
| is given by writing in the numerator the factorials of the different r’s and in 
the denominator the factorials of all the different entries and the factorials of 
_all repeated rows. Thus the coefficient of 


210 
111 is 
111 


413!2! 
201/72! 
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In case two or more functions have the same weight additional equivalent 
partitions are formed by interchange of columns. The reader is referred to the 
above papers for rules for determining the coefficients in the more involved 
cases though the coefficients are presented for all the two way partitions of the 
next section. 


An alternative method of finding the coefficients is that given by C. C. 
Craig (2, 24-25) since it appears that the symbolic formulae used in the present 
paper are essentially his formulae for v’s in terms of \’s. For example his for- 
mula for v4, (2, 22) is given symbolically by the formula for 44 in the next 
section. The only difference revealed is that the subscripts of the \’s are read 
by rows rather than by columns and that they are sometimes interchanged. 
The more precise formulization is needed for the present interpretation although 
it is not needed for Prof. Craig’s purpose. 

A third method utilizes the symbolic multiplication process stated in sec- 
tion 22. Subscripts of the b’s are used to indicate which power sums are col- 
lected. Thus [bo(2) + bu(1)(1))° gives 


babo(4) + baabea(2)(2) + 2[2banbu(3)(1) + baoobou(2)(1)(1)] + 2bubu(2) (2) 
+ 4D; :0b101(2)(1)(1) + Birooboo (1) (1)(1)(1) 


where the underscored terms indicate the products given by [b2(2)]’, 2[be(2)] 
[b1:(1)(1)], and [b(1)(1)]’ respectively. This is represented by 


1 1 4 2 2 4 1 
22 20 21 20 1] 1] 10 
02 01 01 11 10 10 

01 01 01 

01 





r o£ , 
The underscored terms are the only ones remaining when yp; = 0. 

This method is especially useful when a large number of formulae are to be 
computed, as in the next section. 


24. The Partition Representation of Formulae of Total Weight < 8. The 
partition representation of ui(f,) when r S 8 are given in section 21. The 
partition representation of the remaining formulae of total weight < 8, which 
do not contain unit parts, are given below. 


22 1 1 2 
22 20 #11 
02 11 


1 3 
30 12 
02 20 
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4 6 3 
30 21 = 20 
12 21 
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3 15 10 1 15 30 10 5 
51 42 33 50 41 32 23 «40 
02 11 20 O08 12 «21 


60 90 15 30 30 90 

22 13 ~~ 3i 30 12 
11 20 20 20 21 21 
11 11 20 = 02 02 20 


12 16 48 16 
42 33 32 31 
02 ii 12 13 


96 36 48 72 
31 22 30 21 
11 20 12 21 
02 02 02 02 
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6 2 18 
302 212 
030 120 


72 36 
210 111 
111 111 
011 110 


32 3 24 
2111 2200 2011 
0111 0022 0211 


48 48 
2011 
0011 
0200 


24 48 96 16 48 16 32 
2001 2010 2100 0111 1011 1011 O111 
0201 0201 0111 O111 1110 0111 1101 
0020 0011 0011 2000 0101 1100 1010 


1 12 32 12 48 
2000 2000 2000 1100 1100 
0200 0200 0101 1100 0110 
0020 0011 0110 0011 0011 
0002 0011 0011 0011 1001 
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25. The Formulae for the Sample Moments about a Fixed Point in Terms 
of the Moments of the Universe. The partitions of section 21 and section 24 
can be immediately interpreted to give the formulae for the moments of the 
sample function. For example 


un(fs, fe) = bsbonus + (babe + 3baib2 + 6babu)n(n — 1) usue 


/ . . . . . 
and the value of uoi(fs, fe) as given in section 15 can be read by inspection. 
The value of the b’s are to be inserted for any specific function. The coeffi- 

_ S- ‘ ’ i 
cient of us in the expansion of y3(fe) is 
3 2 3 
(bs + 6bebi1 + 8bi1)n(n — 1)(m — 2). 
n— 1 —| . ; 
In case fo = me, bo = ———, and bi: = —~ so that the coefficient is 
n? n 


(n — 1) (n — 2)(n® — 3n” + 9n — 15) 


n> 





as indicated previously by Tchouproff (10, 192) and Church (9, 82). 
The partitions of section 21 give the 8 formulae u,, (x) which Tchouproff 
gave (10, 155). In this case f, = m; and every b is 0 except those having single 
. ] 
subscripts and these equal 7 
The partitions of section 21 give the formulae v,, (7) which were given by 
Tchouproff (10, 186). In this case it is only necessary to take f, = m, and to 
give the b’s the proper values. Tchouproff has arranged his results according 


to decreasing powers of n. As an illustration we derive his result for v4, (~~) = 
! a ‘ 
ui(m,). From section 21 


ui(fs) = banuy + Bdoon(n — 1)u3 
and from Table II 


(n — 1)(n’ —3n + 3) 


by = and be». 





so that 


' 4 
wi(ms) = (1 —-+ = 


n n° 


1 acai — 
= wy t+ : (6u; — 4s) — 2 (15us — 6s) + : (Qu. — 3y,) 


d n 


as indicated by him. 
The partitions of section 24 also give formulae which have appeared before. 
For example the partitions 
1 2 
22 11 
11 
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which symbolize the formula 


uo(fo) = bonus + (b2 + 2bii)n(n — 1)ys 


become 


(n — 1) 


n3 





ys(™me) = [(m — 1)ps + (n® — Qn + 3)y5] 

which was early derived by “Student” (8, 3) and Tchouproff (10, 192). Simi- 
larly the partitions of 222 and 2222 give the formula for and u3(m2) and y4(me) 
which were given by Tchouproff (10, 192-193) and Church (9, 82). 

Sections 21 and 24 can then be used to write the moments about a fixed 
point of a sample function in terms of the moments of the universe. In the 
case of new functions the b’s must first be determined. Formulae involving 
unit columnar partitions are not included. If the formulae were desired in 
terms of moments about a fixed point of the universe, it would be necessary 
to write in addition all possible partitions. See for example the last formula 
of section 23. 


26. The Formulae For Moments of Any Sample Function in Terms of Mo- 
ments of the Universe. The partitions of sections 21 and 24 are also useful in 
writing the formulae for the moments of the sample moments. It is necessary 
to make the usual adjustments in changing from moments about a fixed point 
to moments: 


uo(fr) = ua(f.) — ar (fr) 
wil fr, Sirs) = mir(fr, » fro) sc 110 fer» Sro)Ho1 Sry » Sra): 


The particular two way partitions which are involved in this adjustment are 
immediately recognizable. They are the ones which have an entry which is 
the only entry in the row and in the column in which it is. Thus 3 gives 
220 
002 
one of the terms contributing to u2(fo) 111(fo). In addition its coefficient is the 
same, if sign is not considered, as the coefficient of u2(fo) ui (fo) in the expansion 
of u3(fe) in terms of moments of fe. This has to be true since each is the number 
of ways of forming 220. And so in general the remaining function of n accom- 
002 
panying this adjustment is the product of the coefficient associated with 22 
and that associated with 2. The sign is plus when odd numbers of moments 
are multiplied and minus when even numbers of moments are multiplied. 
Hence 3 contributes —3n°b: to the adjustment to moments and the total 
220 


002 
contribution of 3 to the value of u3(f2) is 3b3[n(n — 1) — n°] = —3b3n. More 
220 
002 
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extensive study leads to the following general method of using the formulae of 
section 24. 
A. Write the coefficient of every two way partition according to section 25. 
B. Block off each single entry by drawing a line through its row and column. 
For example 


6 
220 
602 


The resulting partitions, 22, 2, 2 are called component parts. 
C. Form new partitions by eliminating component parts one at a time, two 
at a time, three at a time, etc. from the original partition in all possible ways. 
D. Form the coefficient of the resulting parts according to the methods of 
section 25. Multiply by (—1)*” where s is the number of resulting parts. 
The values of b will not change. 


E. Multiply in addition by s — 1 when the component parts are all taken 
separately. 
6 
As an example we find the contribution of the partition 2200 to the value 
0020 
0002 


of us(fo). It gives 
6b3[n(n — 1)(n — 2) — 3n°(n — 1) + 2n*Jusuoue = 12ndousud. 


Similarly 1 contributes 
2000 
0200 
0020 
0002 


bs[n — 4nn® + 6n3(n — 1) — 3n‘]us = 3b2(n — 2)uo. 


We use the method in finding the coefficient of y2 in the expansion of u3(mz). 
We find first the coefficient of u3 in the expansion of u3(fe). It is indicated by 
the partitions 


1 6 8 
200 200 110 
020 011 011 
002 011 101 


so that the coefficient of u3 is 

bs[n(n — 1)(n — 2) — 3n?(n — 1) + 2n*] + 6babii [n(n — 1)(n — 2) — n(n — 1)] 
+ 8biin(n — 1)(n — 2) = b3(2n) + 6bebii(— 2n® + 2n) 

+ 8biin(n — 1)(n — 2). 
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° n— 1 —] ; 
When b, = ——— and by = —~ this becomes 
n 





2(n — 1)(n> — 12n + 15) 
~ ~ as 
n® 

previously given by such authors as Tchouproff (10, 194), Church (9, 82), 
Carver (Richardson) (11, 271). 

The general Tchouproff-Church formulae for the third and fourth moments 
of the variance may be written out in this way as may many other moment 
formulae which have not been printed. 


27. The Thiele Moments of the Sample Function in Terms of the Moments 
of the Universe. It is possible also to write the Thiele moments of the sample 
function in terms of the moments of the universe. The technique is very 
similar to that of the previous section. The basis of the transformation is 
now the formula for Thiele moments in terms of moments about a fixed point 
rather than moments in terms of moments about a fixed point. The results 
are the same as those of the last section when a double or a triple product of 
f’s is involved, but they differ with the introduction of a larger number of 
products. The partitions having component parts are broken up into these 
component parts as before but the parts are combined in all possible ways. 
Multipliers are determined as before with the exception that there is a multi- 
plication by (—1)* '(s — 1)! where s is the number of resultant parts. Thus the 
2000 

term 0200 contributes b3[n“ — 4nn® — 3n?(n — 1) + 12n*(n — 1) — 6n'‘]u3 = 
0020 
0002 

—6b3ny3 to the value of \4(fo). 


28. The Moments About a Fixed Point of the Sample Function in Terms 
of the Thiele Moments of the Universe. We return to the problem of section 
25, only we wish to express the results in terms of the Thiele moments of 
the universe. We must use the formulae of section 12. 


r 


1 
‘i z z (Ap,)7! +++ (Ap) 
Pp cee Ps - 


where p; *& 1. 

Thus yu, will contribute to all partitions of r and inversely the contributions 
to a given partition are composed only of these terms which are obtained by 
combining the different elements of the partition. Since the numerical coeff- 
cient in the expansion of yu, is the number of ways in which the r units can 
be collected to form the partition, it follows at once that the complete \ coeffi- 
cient can be obtained by grouping the parts of the partition in all possible 
ways, determining the coefficient of each according to the methods of section 25, 
and adding. In this way the formulae of section 21 can be used to give expan- 
sions in terms of partition moments. For example the representation of u1(fe) 
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1 15 
6 4 
2 


gives at once 
bende + 15[ben + bion(n — 1)JAsd2 + 10[ben + basn(n — 1)]A3 
+ 15[ben + 3bien(n — 1) + Dovon(n — 1)(n — 2)]A3. 


The partitions of section 21 can be made to give the formula y;(l,) which 
were given by Thiele (1, 45-46). For example the formula for yi(f;) is indi- 
cated by 

1 3 
4 2 
2 
so that 


wi(fs) = dandg + 3[bin + been(n — 1))3 


and since 


” (n — 1) = 6n + 6) oat Meee == 
n n 


bs 





(n —1) (nr — 6n + 6)\, 6 (n — 1)r3 
n3 n? 





u3(Is) = ’ 
which agrees with the result as given by him (1, 45). 

The two way partitions of section 24 can be used similarly. This device 
for changing to the )’s is due to the ingenuity of R. A. Fisher who applied it to 
the case where f, = k,. 

As an illustration we write from section 24 the value of u2(fe) in terms of }’s. 
The partition representation 


gives at once 
band, + [bon + dan(n — 1)]A3 + 2[bon + diun(n — 1)]d3 


which agrees with the result of section 12. The other illustrations of that 
section may be written out similarly. 

As a final illustration of this technique we find the coefficient of \j in the 
expansion of ya(fs, fe). The partitions are 


2 9 18 
301 220 211 
031 112 121 
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and the coefficient is 
2[bsben + b3bun(n — 1)] + O[bsben + dirben(n — 1)] 

+ 18[b3bon + beibun(n — 1)] + 6[b3ben + dsbeiben(n — 1)]. 
If the b’s are inserted to form the k’s, the first and last terms become 0 and the 
2in — 45 


others give - 


n(n — 1)? 
(3, 208). 





This agrees with the value as given by R. A. Fisher 











29. The Moments of the Sample Function in Terms of the Thiele Moments 
of the Universe. The partition representations of section 21 and section 24 
can be used similarly to write formulae for the moments of the sample function 
in terms of the Thiele moments of the universe. It is only necessary to use the 
general plan of section 26, but to write the coefficient of every resulting parti- 
tion according to the method of section 28. For example the partition 


20 






gives the coefficient 





bs[n + 4n® + 3n® + 6n® + n®] — 4b3[n? + 3n°(n — 1) + n(n — 1)(n — 2)] 
+ 6b3[n*> + n°(n — 1)] — 3bon* = ba[n* — 





4n* + 6n* — 3n'‘] = 0. 

















30. The Thiele Moments of the Sample Function in Terms of the Thiele 
Moments of the Universe. The partition representations of section 21 and 
section 24 can also be interpreted to give the Thiele moments of the sample 
function in terms of the Thiele moments of the universe. The scheme is 
similar to that of section 29 except that the formulae for changing to Thiele 
2000 
moments are used as in section 27. For example the partition 0200 has now 
0020 
0002 


associated with it 
ba[n + 4n® + 3n™ + 6n™ + n™] — 4b3[n? + 3n7(n — 1) + n(n — 1)(n — 2)] 
— 3b3n’(n — 1)” + 1263 [n® + n3 (n — 1)] — 6ban* = 0. 


The application of this method enables one to write the formulae of section 13 
(and others which they typify) with relative ease. It is now possible to com- 
plete the task left unfinished in section 15. We do not take the space necessary 
to write all the terms of da(f3, fe) since the lengthy expression can be obtained 
quite readily from the representation of section 24. One term, say the coeffi- 
cient of AsA2, Is represented by 
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1 9 12 6 
330 222 321 
002 110 011 


and gives 
Q[bsbon + boibon(n — 1)] + 12[b3ben + dbsbabun(n — 1)] 
+ 6[b3bon + dsbaiben(n — 1)] 


—|} 


n(n — 1) T his 


which becomes ae when b3 = bo = 1 and by = by = 
n(n — 1) n 
agrees with the result given by R. A. Fisher (3, 209). 

For simplicity of form it is logical to use this formulization of results, Thiele 
moments in terms of Thiele moments, and it has been used by Thiele (1), 
Craig (2), Fisher (3) and Georgescu (4). They however have used different 
sample moment functions. Thiele and Georgescu used the Thiele moments 
of the sample, Craig and Georgescu the moments while Fisher introduced the 
k function. 

The present discussion deals with the corresponding partition moments of 
any rational integral isobaric moment function of the sample. The results 
indicated here give many of the results of the previous authors as special cases. 
For example the symbolic formula 44 of section 24 gives the m’de(us) of Thiele 
(1, 45), the Soe(ve, vs) of Craig (2, 57), the «(44) of R. A. Fisher (8, 210) as 
special cases when the formula 44 is given the interpretation of this section. 

Some may prefer the Craig attack (2, 21-35) to the partition method. It 
should be noted that the formulae of sections 21 and 24 can be used in place 
of part of the Craig method. Thus his formulae (2, 22) 


Vso = Aso + 28 Agsorzo + 56 Asorzo + ete. 
vag = Ags + (12 Asdo2 + 16 AssA11) + etc. 


are immediately obtainable from the symbolic formulae by writing )’s in place 
of b’s and by using row, rather than column, subscripts. It is then necessary 
to compute the values of \j,;. --. as given by him (2, 16-17, 40) and to insert 
in his expansions of Sy(vm, vn) in terms of v’s. For example 


Su(v3, v2) = . [v50 + (n = 1) v32 = Nv30V02] (2, 32) 


and from the symbolic formulae of sections 21 and 24 
¥50 = Aso + 1LOA30Az0 
= Aze + AgoAo2 + SAi2A2%0 + 6A21A11 
i ts 
20 
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so that 





] 
Si(vs , v2) = . [Aso + (n sie 1)Aze + QAz0A20 + (n — 1)(6daA11 + 3AcrA20) | (2, 30) 


which agrees with that given by Prof. Craig (aside from an obvious typographical 
error). The insertion of the values of gives the value as indicated by 
Au(ms3 , M2) of section 13 and by the first method of the present section. 





31. Special Rules for the Determination of the Coefficients in the Case of 
the Fisher and Georgescu Analyses. R. A. Fisher (3) gave a number of simple 
rules which assist greatly in the determination of the coefficients accompanying 
the partitions. Georgescu (4) also introduced special rules for the evaluation 
of the coefficients of the different partitions he used. It is not to be expected 
that all these rules are applicable in the more general case under present con- 
sideration, but the vanishing of such coefficients as that of 2000 leads one to 
0200 
0020 
0002 
suspect that there might be some rules which are applicable to this general 
case. A sensible method of procedure is to examine the rules of Fisher and 
Georgescu and determine if they hold in the more general analysis. The special 
rules of R. A. Fisher might be given somewhat as follows. 
A. If a partition has a column with a single entry, that column may be 
eliminated and the factor n~' introduced. 
B. Any partition having a row with a single entry may be neglected. 
C. “We may exclude any partition in which any set of rows is connected 
to its complementary set by a single column only.” 
D. In determining the algebraic coefficient of a partition the “‘pattern’’ is 
sufficient and precise entries are not needed. Thus the partitions 21 and 35, 
11 42 
although they have different numerical factors, have associated with them the 
same function of n. This value is indicated by the pattern zx which has asso- 
rx 
















: iis ; | , . 
ciated with it the function scene As a result of this property Fisher was able 
n —_— 


to provide a table (3, 223-226) of useful patterns which is of great assistance 
in writing the value of the coefficients. 

E. Formulae of moments of k functions involving k; can be derived from 
corresponding formulae not involving ki. “The effect upon the corresponding 
formula of adding a new unit part to the partition is (1) to modify every 
term in the formula by increasing the suffix of one of its x functions by unity 
in every possible way, and (2) to divide the whole by n.”’ (3, 206). 

Two of the important Georgescu rules may be stated. 

A’. The numerator function (aside from numerical coefficient) is not altered 
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° : 1 . 3- 
if columns are changed to rows and vice versa. Thus the coefficient of s2 in 


6N(N — 1) 4N(N — 1) 


S(3°) = WW +)? and the coefficient of s3 in S(2°) is ——— . Georgescu 
has replaced n by N + 1. 

B’. All partitions which can be broken up into component parts have coeffi- 
cients of 0. This is extended to include all partitions which have as component 
parts other partitions. Thus 


a + iP 


2100 
1100 
0012 
0034 


has a coefficient 0 as does the equivalent 


2010 
1010 
0102 
0304 


32. Special Rules for the Determination of the Coefficients in the More 
General Case. In the more general case we have 

A. If a partition has a single column with a single entry, c, that column 
may be eliminated and the value b, inserted as a multiplier. This is imme- 
diately evident since the contribution of that column to each term in the 
expansion is b, times its value if the column were eliminated. 

B. The coefficient of any partition having an entry which is the only entry 
in its row and column, is 0. 

This rule, which saves considerable labor in that it makes unnecessary the 
computation of the coefficients of many of the partitions of section 24, is estab- 
lished in this way. Without loss of generality the partition may be repre- 
sented by 


Cu 
C21 


Tuti,vtl = C31 


0 


Cu+i1, v+l 


and 7, ,» may represent the partition containing the first « rows and the first v 
columns. We determine the coefficient of 7.41,,41 in terms of the coefficient 
of r.,,». Consider first any grouping of the uw rows of 7, ,, into w rows. There 
will be w corresponding groupings of 7,41,.41 in which the last row is added, in 
turn, to each of the w rows and another w + 1 rowed term in which it is not 
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added. In each of the first w cases the coefficient by rule A is multiplied by 
Beus1,e4:- In the case of the w + 1 rowed partition the coefficient is multi- 
plied by b.,,,,,..;, and n“” is replaced by n‘“*”. A final adjustment takes 
care of the transition from the moment about a fixed point of the sample 
function to the Thiele moment of the sample function. This adjustment de- 
mands the multiplication of the coefficient of 7, ,. by b.,.,,,4, mand the sub- 
traction from the sum of the other terms. If B, is the coefficient of the w 
rowed form, it follows at once that the corresponding coefficient is 
















(w) »+1) (w) 
Budex, 41 (wn +n” — nn] = 0. 


This holds for the expansion of any term of z,,, and hence the coefficient of 
Tu+1,r41 is 0. Of course the argument holds if the partition has more than 2 
component parts. 

It thus appears that this rule holds not only for k, and m, as Fisher and 
Georgescu have noted, but for f,. 

C. The coefficient of any partition which can be broken into component 
parts is 0. In this sense a component part is any group of rows or columns 
which have no entry in common with any other group of rows or columns. 
It corresponds in matrix language to a matrix which results when one matrix 
is zero bordered by another matrix although rows and columns may thereafter 
be interchanged. 

The proof of this more general case follows the general line of the simpler 
case although the reasoning is more complicated. For example the coefficient of 








Cn C++: Cw O 0 


0 










C9 + > 






0 


Cu2 °° Cu 0 









0 ee 0 Cu41, 041 Cut , v42 






0 0 T= 0 Cu+2 »v+l1 Cu+2 , v42 









is 0 since any w rowed term of the 7, ,, contributes 
(w) (w+1) (w) 
Bude, + > vtitCut+2» etl bess svt2tCut+2 » v+2 [wn + n — mn ] 
(w) (w+1) (w+2) 
ee eee ee we +n 
(w) 
—n(n—1)n”]=0. 


Other special rules of Fisher and Georgescu do not hold in the general case. 
Thus Fisher rule B is not generally true since the partitions 


12 
30 








22 
20 


and 
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have respective algebraic coefficients of bsbon + bsiben(n — 1) and 
babon “+ booben(n = 1) 


and these are not in general equal to 0. 

The Fisher rule C is replaced by the somewhat less general C of the present 
section. 

The Fisher rule D is not applicable in the general case. The Fisher rule D 
is applicable in all cases in which the value of the 6,7: ... ,%. is completely deter- 
mined by the number of parts for in this case the particular value of each 
part is not pertinent. We may say then that the Fisher rule D is applicable 
to all cases in which 6,71,%2 ... ,%» is a function of p, n where p is the number 

ep 
of parts. This condition is satisfied by 6,71... ,%° = —— and the 
coefficients are worked out for it in Fisher’s paper. The same method is 
applicable to other functions satisfying the general condition although the 
values of the coefficients will of course vary with the definition of b. 

The Fisher rule E is not applicable to the general case. Its validity, from 
an algebraic standpoint, depends upon the Fisher property B which is not 
generally applicable. The Fisher rule E as applied to the more general case 
gives correct terms but it does not give all the terms. For example the Fisher 
rule E applied to A2(ke) gives 


r Qr3 
Ne(ke) = : n | 


4dsdo 
n(n — 1) 


r 
Aai(ke, ki) = = + 


The application of a corresponding rule to 
Ao(fo) = bendy + 2[bin + bun(n — 1)]A3 
would give 
or( fo, fi) = babinds + 4[bobin + Biibin(n — 1)]Asre 
while the correct result is indicated by 


1 4 
221 210 
011 


and is 
hai(fofi) = bebinds + 4[bobin + bobubsn(n—1)]As\2 + Q[bsbin + bbin(n — 1)]AsA2 
+ 4[b3bin + biibin(n — 1)]Asdz. 
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The difference is due to the vanishing of the two middle terms in the case of 
the k functions. 

The rule B’, which Georgescu found most useful in computing and checking 
his formulae, is not generally true. It is not even true in the case of the k 
function, as can be discovered by using it on the list given by R. A. Fisher 
(3,210). It is interesting to note that the Georgescu method, while not being 
able to utilize many of the special rules of the Fisher method, does use this rule 
which is not in general adaptable to the Fisher method. 





33. Special Rules in the Case of the h’ Functions. Special rules can be 
worked out for other sample functions. As an illustration we examine the 
: : . : : 1 
function h. which was defined in section 19. It is recalled that by = os and 

n 
that b,7:... ,%» = 0 for all other cases. It follows at once that 
A. Any partition having any entry other than unity (or zero) may be 
neglected. 


B. The value of by is i 


np) 





‘ * ‘ , , ie a 
As an illustration we write the value Az(h3, he). From the partitions of 
section 24 we select 


36 36 
111 110 
111 and 110 
110 10] 


O11 





as being the only partitions making a contribution. The result of section 19 
follows at once. 


34. The Case of a Normal Universe. A normal universe is characterized by 
the relationship that 4, = 0 when r > 2. It follows that it is only necessary 
to compute the coefficients of those partitions giving powers of d2. 

Wishart (5) (7) has developed the partition analysis of the k function in 
the case of a normal parent while Georgescu has studied the corresponding 
m function. It is not the purpose of this section to make extensive study of 
the case of the normal parent but simply to indicate that the results of section 24 
are immediately applicable. As an illustration we write the values of Ai(f2), 
A2(fo), As(fe) and A4(fo) in the case of a normal universe. The terms are given 
successively, by 


1 2 8 48 
2 1] 110 1100 
11 O11 0110 

101 0011 


1001 
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and hence 
hi(fe) = bende 
ho(fo) = 2[bsn + biin(n — 1)]A5 
Aa(fe) = 8[b2n + 3bebiin(n — 1) + bin(n — 1)(n — 2))r2 
ha(fe) = 48[bon + 6b3biin(n — 1) + bin(n — 1) + 4bobin(n — 1)(n — 2) 
+ 2biin(n — 1)(n — 2) + diin(n — 1)(m — 2)(n — 3)]A3. 


It is only necessary to substitute the b’s to obtain the results for different values 
of f. This is done in Table IV. 


TABLE IV 
The first four Thiele moments of fe for various sample functions in the case of a 
normal universe 


Sample 


func- Ai (fo) 2 As( fo) | dal fe) 
tion 





(n-1), 


nm 


Ae 


(n —1 ) Ne 
n 





nee we | 





8(n — 1) .3 
on re 
__ 82 
(n — 1)? 
8(n — 1) 


n§ 
ni 
8d2 

(n — 1)? 








n(n — 1)? 


8(n — 2)r3 | 





48(n — 1) 2 
nt 
482 _ 
(n — 1) 
48(n—1)d2 
n4 
483 
nt 


_ 482 


(n — 1)3 
48(n” — 3n + 3)A3 


n(n — 1)8 





One surmises that the general value of 


d, (fe) is 2° "(r — 1)! A2B: 11000 --- 
01100 - 


00000 --. 
10000 --.- 
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where B represents the b coefficient of the r rowed partition. This induction 
appears consistent with the fact that 
o'r! — 
Args(ke) = ——— 
(ha) = 
as shown by John Wishart (7). The whole subject of the Thiele moments of 


the general function in the case of a normal universe would make an interesting 
subject of investigation. 






















35. Summary and Conclusion. The contributions of this paper include 
1. The definitions of specific moment functions in terms of power sums. 
2. The use of indeterminate multipliers in representing a general isobaric 
moment function. 


3. The finding of the expected value of products of these functions by alge- 
braic methods. 

4. The use of tables in writing these expected values in terms of moments 
(or of moments about a fixed point) of the universe. 

5. The finding of the expected values of specific moment functions by sub- 
stitution. 

6. Means of establishing the expansion of new moment functions which are 
defined by their expected values. 

7. The introduction of the sample function of weight 7 whose expected 
value is p,. 

8. The introduction of the sample function of weight r whose expected 
value is p,’. 

9. The two way partition formulae of weight < 8 which do not involve 
unit parts. 

The use of these partition formulae in writing: 

10. The moments about a fixed point of f, in terms of moments. 

11. The moments of f, in terms of moments. 

12. The Thiele moments of f, in terms of moments. 

13. The moments about a fixed point of f, in terms of Thiele moments. 

14. The moments of f, in terms of Thiele moments. 

15. The Thiele moments of f, in terms of Thiele moments. 

16. Special rules in the case of Thiele moments. 

17. The applicability of these results to a given sample moment function 
and hence the derivation of varied results, of such authors as Thiele, Tchouproff, 
Church, Fisher, Craig, and Georgescu, from the same partition formulae 

18. The simplicity of the formulae when hi. is used as the sample function. 

19. The application of the synthetic formulae to the Craig method. 

20. The applicability of the theory to a normal universe. 

The introduction of such general procedure opens up a wide field for future 
study. It is impossible in a single paper dealing with so broad a subject to do 
more than to outline the general scheme by which two way partitions can be 
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used as a central formulization of the various formulae for moments of moments. 
More detailed proofs and more extensive analysis of the more important of the 
special cases will undoubtedly be supplied by later writers. 

In later papers the author will show how the partition representation can 
be used in the case of multivariate distributions and how it can also be used, 
in connection with the sampling polynomials introduced by H. C. Carver (11), 
to represent the more complex formulae obtained in the case of finite sampling. 

It is obvious that the author is indebted to the classical moment studies of 
Fisher and Craig. He also wishes to acknowledge his indebtedness to Prof. 
Craig and to Prof. Carver who have read the manuscript and have made 
valuable suggestions. 


Tue UNIVERSITY OF MICHIGAN. 


REFERENCES 


(1) Turete, T. N. ‘‘Theory of observations.’”? London 1903. C. and E. Layton. Re- 
printed in ‘‘Annals of Mathematical Statistics,’ 2 (1931), pp. 165-306. 

(2) Craic,C.C. ‘An application of Thiele’s semi-invariants to the sampling problem.’’ 
Metron, 7 (1928-29), pp. 3-74. 

(3) FisHer, R. A. ‘‘Moments and product moments of sampling distributions.”’ 
Proceedings of the London Mathematical Society. Series 2, vol. 30 (1930), 
pp. 199-238. 

(4) Sr. Georcescu, N. ‘Further contributions to the sampling problem.’’ Biometrika, 
24 (1932), pp. 65-107. 

(5) WisHart, JoHN. “The derivation of certain high order sampling product moments 
from a normal population.’”’ Biometrika, 22 (1930), pp. 224-238. 

(6) Soper, H. E. ‘Sampling moments of moments of samples of n units each drawn 
from an unchanging sampled population, from the point of view of semi- 
invariants.’”’ Journal of Royal Statistical Society, 93 (1930), pp. 104-114. 

(7) WisHart, JoHN. ‘‘A comparison of the semi-invariants of the distributions of 
moment and semi-invariant estimates in samples from an infinite population.”’ 
Biometrika, 25 (1932), pp. 52-60. 

(8) ‘‘Student.’”’ ‘The probable error of the mean.’’ Biometrika, 6 (1908), pp. 1-25. 

(9) CHurcu, A. E. R. ‘On the moments of the distribution of squared standard devia- 
tions for samples of » drawn from an indefinitely large population.’’ Bio- 
metrika, 17 (1925), pp. 79-83. 

(10) TcHouprorr, A. A. ‘‘On the mathematical expectation of the moments of frequency 
distributions.’’ Biometrika, 12 (1918-19), pp. 140-169, 185-210. 

(11) Carver, H. C. ‘‘Fundamentals of sampling.’’ Editorial. ‘‘Annals of Mathe- 
matical Statistics,’’ 1 (1930), pp. 101-121; 260-274. 

(12) Frisco, Ragnar, ‘‘Sur les semi-invariants et moments employes dans |’etude des 
distributions statistiques.’’ Skrifter utgitt av det Norske Videnskaps-Akademi 
i Oslo. II Hist.-Filos Klasse (1926), no. 3. 













NOTES 


A COEFFICIENT OF CORRELATION BETWEEN SCHOLARSHIP 
AND SALARIES 


INTRODUCTION 


Some might doubt that it is correct to apply a coefficient of correlation to 
show the relationship between scholarship and salaries. This coefficient can 
be trusted to give at least a rough approximation, which is all that is necessary 
in the inexact science of vocation. It is fictitious accuracy to be too finical 
in the application of formulas. Therefore, a coefficient of correlation between 
scholarship and salaries is a valuable part of human knowledge. 

Would it be worth while to find this coefficient if it is based upon the experi- 
ence of the American Telegraph and Telephone Company? Since the employ- 
ment practices of this company are not representative of the employment 
practices of business at large, one might doubt the validity of drawing general 
conclusions from such specialized data. The coefficient for business at large 
is probably less than the coefficient for the Bell System; the value of this knowl- 
edge is enhanced if we know the latter coefficient. Since this company is very 
large, a coefficient between scholarship and salaries would be valuable, even if 
this coefficient applies only to the Bell System and to other companies having 
approximately the same employment practices. 

An article’ by Mr. Walter 8. Gifford, President of the Bell System, contains 
a discussion of some of the relationships between scholarship and salaries. 
President Gifford, however, did not determine in the case of the Bell System a 
coefficient of correlation between scholarship and salaries. 

The purpose of this article is not a new contribution to statistical method, 
but is an application of the method’ of finding the coefficient of correlation 
when the two variables have not been quantitatively measured. This method 
will be applied to the chart on page 672 of President Gifford’s article, in order 
to determine for the Bell System the coefficient of correlation between scholar- 
ship and salaries. 




























FINDING THE COEFFICIENT OF CORRELATION 


An explanation of the chart. It is based on the experience of 2,144 Bell 
System employees over five years out of college. First, assume these employees 





1 It is entitled ‘“Does Business Want Scholars?” and was printed in the May 1928 issue 
of Harper’s Magazine. 
2 It can be found in Elderton’s ‘“‘Frequency Curves and Correlation.”’ 
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are grouped according to their grades in college. In the high scholarship group 
put those who graduated in the highest third of their classes. The middle 
and low scholarship groups are formed in like manner. Secondly, suppose the 
same employees are divided into three equal groups according to their salaries. 
Then, the salary of any one of the employees would be high, middle, or low. 

Assume a hypothetical group of 300 employees who are college graduates. 
Suppose that the scholarship of 100 of them was high, that the scholarship of 
100 of them was middle, and that the scholarship of the others was low. Also 
assume that the salary experience of these 300 employees is the same as that 
of the 2,144 employees of the Bell System. 

The 300 employees can be grouped according to the following table. 


TABLE NO. 1 





| Scholarship 
_o eS naan Totals 
Middle | 





39 





This table can be combined as follows. 


TABLE NO. 2 


Scholarship 





Salary 


Low & Middle | High 


| c | d 


Middle & Low a b 





Then, c = 46, a = 154, d = 48, and b = 52. Assume N = 300. 

Assume 2 is a function of grades received in college. Suppose y is a function 
of salaries received. Assume that the frequencies x and y both follow the 
normal curve of error whose standard deviation is equal to one. Also assume 
that the average of x and the average of y are both equal to zero. It is a 
matter of common knowledge that salaries are not arranged in a symmetrical 
fashion; y is not a linear function of salaries. 

In the formulas which follow, r is the symbol for the coefficient of correlation. 
These formulas are applied to Table No. 2. We have 

1 [ ot? ae (a+c) —(b+d) 

0 


—— dx = ——__-___ = .167, andh = .4316. 
/ 2x 2N 





k 
= I et dy = aR = .187, andk = .4874. 
Tr J9 


Then, 


1 —1p2 y 1 
H = ——e”™ = .3635, and K = 
V 21 V 2r 


All the quantities except r in the following approximate equation are known: 


ad — be 
N?HK 


e ** — 3543. 


2 3 
=r+5hk+ (hr — 1) -1) 


4 5 
i 2 oe r 4 _ gp? 4 gz2 
+ 57 Mh’ — 3)k(R — 3) + Fo (ht — 6h? + 3) (RY — 6B? + 8). 


Therefore, 
0261r° + .0681r* + .1034r*° + .1052r + r — .4314 = 0. 


Then, r is approximately equal to .4051. Consequently, for practical purposes 
we can assume that r = .4. 
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NOTE ON THE DERIVATION OF THE MULTIPLE CORRELATION 
COEFFICIENT 


Consider N observed values of each of n variables. These n-N values may 
be tabulated in a double-entry table as follows: 


Xu Xi Xis--- Xin 
Xo X22 X23 2 X N 


Xu X ne X n3 ee Xan 


where Xx is the k** value of the 7*® variable. 
Using the i” variable as the dependent variable, the general linear relation- 
ship between the n variables may be expressed by 


i= H+ mt --- + Gatat Gntig +---+ drt, (1) 
where 

ia; is the general parameter which is to be determined empirically; 

a; = X; — M;; 


M; is the arithmetic mean of the j'* variable. 
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By the method of least squares, the constants of (1) must satisfy the normal 
equations: 
(E23) ir + (24x42) de + --- + (2x41) as-1 
+ (Satis) Ginn +--+: + (2x22), = 
(242%) 1 + (S25)ie + > + (Laer ;-1) Ais 


+ (Sxer ins) Giga +--+ + (Lretn) in 


(224-171), + (VAi1T2) ie + --- + (LAiuT,)a. = 


(Sais) Gy + (Sa ipite) de +--+ + (LHiFn) in = 


(Stat), + (VInT2)2 + --- + (Tr%)idn 
where 


. 
(Sozz;) = 2 (Xun — M,) (Xn — Mj). 
(L2,2;) = Nri;0;0;, 
Sri) = Noi = Nruoww: (2) 
where 
r;; is the Pearsonian coefficient of correlation between the 7*” and j** variables, 
o;, the standard deviation of the 7** variable. 
Substituting the right members of (2) in the normal equations, we obtain 


the system: 


n 
> T1K010% ian = 
k=1 


n 


> Tox020k ik: 
| 


n 


7 Ti-1,k Gi-10% x = O 
k=1 


n 
Zz Tina, k C1410; 1, = O 
k=} 


n 
y » TrkOnd At = 0 
mi 
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where 


(4) 


A,; be the first minor of the element r;jo;0; in A, A be A with the z** and kth 
columns interchanged, and ;,A;; be the first minor of the element in the 7‘ 
column and 7** row of 4A. 

Solving (3) for ,a, by Cramer’s rule, we find 


But it can easily be proved that 


uA = (- 1)" Aa; 


hence 


et Ax 
iat = (—1) , ms 


Using cofactors of A instead of minors, we have 


L)"Dx _ _ Du 


ee (— 


iat = (- 


Di - Di 


Without writing the determinant out in full, we notice that the o’s can be 
factored out. Hence 


9 9 


22 2 2 2 2 2- 
0102 -** Ok-19kOk+1 °° * TO{-19j;0j41°°* On Kix 


e@ 9 9 


9 2 r 
0109 +++ O;-10i41 «++ On Ki 


where 
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Using these derived values for the coefficients, we may write (1) in the sym- 
metric form: 


Pix. ~ mw) 4% “2 (X_— Ma) +--+ + A(X, - M,) = 0, 
1 


o 


be wa ™ 


gs 
For a multiple correlation coefficient, we use the formula 


N i—1 n 2 
ie E - (= A Lnj + a aus) 


R? ai _ im! k=1 k=i+1 


2 
N o; 





which measures the amount of observed dispersion from the regression plane 
in which X; is the dependent variable. 
Substituting the values for the a’s, we find 


(Fae. 5 — +> +i | 


on 


72 
it N 


Squaring the bracket expression and using (2), we obtain 


, , N 
R? KiKi 7 LejL1j 
v 


i=1 
Nox0; 


= = KaKurs| 


k=1 (=1 


1 [Elke] 


i=1 


The second sum is the sum of the products of the elements in the kt® row 
by the cofactors of the elements in the 7** row. This sum is necessarily zero 
unless k = 7; but if k = 7, this sum is equal to K. 
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NOTES 
NOTE ON NUMERICAL EVALUATION OF DOUBLE SERIES’ 


1. The Euler-Maclaurin summation formula has been extended to two 
variables by Dr. Sheppard,” and Mr. Irwin,’ to determine cubature formulas. 
A more complicated two-dimensional form was given by Baten‘ involving 
product polynomials, for which a remainder term was also calculated. The 
purpose of this note is to apply the simpler formula to the numerical evaluation 
of double series of positive terms. The method may be extended to multiple 
series of order p > 2. If the double series converges one may sum by rows 
(or columns), using the ordinary sum formula twice. The method is to take 
out a rectangular block of mn terms and then apply the formula to the remaining 
terms. By taking m and n sufficiently large one may cause the series resulting 
from the formula to converge sufficiently rapidly to obtain the sum to the 
desired number of decimal places. For practical work the error may be es- 
timated because of the asymptotic character of the series involved in the Euler- 
Maclaurin formula. 
Write this in the form 


| © f(z) = [ het = yo 279 , Fe -P'e _ 

] ‘ . ‘ - 9 9 ) 

ry (a) ae sf (s) f N(a) - f a) a “ (_ 1B f° (a) —f* (s) 
30240 1209600 (27)! 














If s— one has accordingly in the ordinary case of convergence 


7 ” / tr V 
(2) dX f(z) = [ f(x)dx + $f(a) — +. j oa fe ie i 








co 


Now define v(x) = >» u(x, y) = [ u(x, y) dy +°3u(z, b) — aS + 
y=b b 
u,a(x, b) _ 


son and wly) = >» u(z,y) = u(x, y)dx + Zula, y) — ue(a, y) y) 4 
720 z=a a 2 


me = a—1 b-1 ~% b—1 
tust(a, ¥) _ --+, then D> Do ulz,y) = DO DY alz,y) + Dd v(x) + D wly) 
720 ie 2 2 2 ~ 





sa ” ; 1. a v’(1) v’"(1) - v’(1) - a—1 b—1 
= / v(x)dx + 4v(1) 2 + 720 30240 + + a a u(z, y) 


(3) 


b , uy? at 
(y)dy — 4w(b) + aw(1) + 20) = wD _ wit) = wD) 
+f w(y)dy — 3w(b) + 3w(1) + 7s - ex's 


1 Presented to the Society, Nov. 30, 1934. 

2 W. F. Sheppard, ‘‘Some Quadrature Formulae,’’ Proc. London Math. Soc., Vol. xxxii, 
1900. 

3 J. O. Irwin, ‘‘Tracts for Computers,’’ No. X, Cambridge Univ. Press, 1923, On Quad- 
rature and Cubature. 

4W. D. Baten, ‘‘A Remainder for the Euler-Maclaurin summation formula in two 
independent variables,’’ Amer. Journal of Math., Vol. 54, 1932, pp. 265-275. 
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a—1 b—1 


aw a—1 
Instead of this one may use a. u(x, y) + i wy) + > v(x). The scheme 
z=1 y=1 y=1 z=1 


of the double series may be illustrated by a sketch of a quadrant of the zy-plane 
in which the point (x, y) represents the term u(z, y). 
Evidently by taking a combination of results from (3) one may evaluate 


q t 
quite readily such finite sums as >> > u(z, y) where q and ¢ are large. 


z=p y=? 

As an illustration of (3) consider >> >> (x” + ay’). Here one needs to 
evaluate the integral of the summand. The transformations z = ay tan 6 
and y = 1/t lead to a form which may be integrated by parts. The more 
complicated form >, >. (ax” + 2bry + cy’) for the case in which s > 3/2, 
a > 1, might be handled by using x = 1/t and approximate integration by 
Simpson’s rule. 

Take as a second example > d“(e@+ y) ”,p > 2. The case of p = 4 was 
carried out by taking a = b = 10 in (8) and carrying the computation to twelve 
decimals. The series involved converge rapidly and a result was obtained 
which differed by 2 in the 12th place from the true value 0.119 733 669 448°. 


By summing diagonally one may convert this to the simple series > 2(z + 1) 
1 


or >, (s—1)s* = D> (s *— ss“). The method of summation diagonally may 
2 ] 


be extended to > >> (x + ay) ’, p > 2,a > 0, by the applications of the 
Euler-Maclaurin sum formulas (1), (2) in succession after a triangular array 
of terms have been omitted. 


The form >> >> 2” y “ can be written as the product of the single series 


(2.2 "XZ 9"). 


2. Another method of numerical evaluation is the analog of that used for 
single series by the author.° Instead of rectangles one has right prisms of 
square or rectangular cross-section. Instead of shifting the rectangles one unit 
to the right to determine upper and lower bounds the prisms are shifted diago- 
nally so that they go effectively one unit in each variable. In the case of a square 
base each prism is moved along the 45° line one diagonal unit length. For 
the lower bound instead of trapezoids one uses truncated prisms. For example, 
the prism of height umn is cut by two planes, one determined by the upper 
vertices Umn, Um, nti, Umsi,n and the other by the upper vertices Umsi,n, Um,n+t) 
Um+i, n41 Of the truncated prism. The surface z = u(m, n) passes through 
all the upper corners of the truncated prisms. Each prism is composed of 
two truncated triangular prisms. Now the volume of such a triangular prism is 
the arithmetic mean of its vertical edges multiplied by the area of its base. 


5**& New Method for Finding the Numerical Sum of an Infinite Series,’’ Amer. Math. 
Monthly, vol. XL, No. 9, Nov., 1933, pp. 537-542. 
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Hence the difference in volume between the truncated rectangular prism men- 
tioned above and the prism of uniform hight z = u,,, can be shown to be 
(4) 


(Stimn = Cesiiwt Bast, 0 —_ 2m, nii)/6. 





Let us consider series whose corresponding surfaces do not rise above these 
truncated prisms. This sort of truncated prism differs less from the volume 
under the surface than the one formed by the diagonal joining the other pair 
of upper vertices and planes through it for upper faces. The lower bound for 
the remainder is the volume under the surface extending to infinity in the 
m and n directions plus the sum of these differences. Accordingly one deter- 
mines as the lower bound for the remainder R»_;, »-1 after summing a rec- 


m—1 n—1 


tangular array » : u;,; the form 


7=1 j=1 


m 


2 oo 
es Ds ~ Se 4 Od ees + 92, ees + 32 


- i=1 j7=1 t=1 
(5) 


n 2 a i) m 
4+ » Umi + / | Um,ndmdn + I i Um,ndmdn < Ry—1,n-1. 
p=1 1 m n 1 


The upper bound may likewise be given as follows: 


(6) | << S + T + I Um,n dmdn ail K| 7 Um-1, j + p vin | 
m—1 d 


n—l j=n—1 =m 


tol 


where 


x n—l1 2° 060ml 
(7) S=) Lui, T=), Duis, 
i=m j=1 j=n i=1 


7 2 2 30 20 co) 3 
Um.rndmdn — Um.ndmdn — > Um—1,j — 7. Nitcont 
(8) k sie n—l m—1 _—, _ Jn mm ae EL __t=m_ Mad 


Um—1,n—1 + Um,n—-1 


An alternate definition of k is 


n m 
(9) ‘” = | | [ Um.» dmdn va vn | = a Pe rv Hic sa 
n—l Jm—1 


An illustration is afforded by Lo Z. (m + 1)°% for which k = .45614, 


n=1 m=1 
k’ = .44536 when m = n = 10 in (8), (9). In this case (5) gave an error of 
—14 X 10° and (6) an error of 10°. 
S and T may be evaluated by the method published in the Monthly.° 
One must assume that k increases with m and n._ It is evident that for this 
























6 Loc. cit. 
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method and for the one in the Monthly differentiability is not required but 
only integrability, conditions less restrictive than those required by the Euler- 
Maclaurin summation formulas. It is also clear that the method may be 
extended to multiple series of positive terms of multiplicity greater than two. 
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REPORT OF THE ANNUAL MEETING OF THE INSTITUTE OF 
MATHEMATICAL STATISTICS 


The meeting of the Institute of Mathematical Statistics for 1936 was held in 
Chicago on December 28-30 in connection with the meetings of the American 
Statistical Association and the Econometric Society. 

In addition to the sessions at which voluntary papers were read, a session with 
invited papers was held on the morning of December 30. At the invitation of 
the Program Committee, Professor P. R. Rider presented a paper on “‘Recent 
Advances in Mathematical Statistics: Factorial Design’’ and Professor Harold 
Hotelling spoke on ‘‘The Analysis of Sets of Correlated Variates.”’ 

Professor C. C. Craig of the University of Michigan and Professor A. R. Cra- 
thorne of the University of Illinois constituted the Program Committee. 

At the business meeting of the Institute, the following officers were elected 
for the year 1937: President, Dr. W. A. Shewhart; Vice-Presidents, Professors 
P. R. Rider and B. H. Camp; Secretary-Treasurer, Professor A. T. Craig. 

The Institute voted that it would presumably hold its 1937 meeting with the 
American Mathematical Society. 

ALLEN T. CraliG, 
Secretary. 


NOTICE TO SUBSCRIBERS 


Plans are under way to include in the ANNALS a new section, entitled ‘“Numer- 
ical Illustrations of Statistical Methodology.” This new section will be a 
regular feature of the ANNALS, and will deal with the application of statistical 
technique and theory to the solution of problems in various fields. It is hoped 
that this new section will be of considerable value to those who are primarily 
interested in numerical applications of the more recent theoretical developments 
in mathematical statistics. 

The Editor will welcome contributions to this new section of the ANNALS. 
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